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SUMMXy 


Ourlng  the  oouree  of  this  oontreot  an  extensive  series  of  propwee 
have  been  prepared  for  an  IBM  7070  data  prpoeaelng  nadilne*  Ihla 
aet  of  prograsu  Inoludes  oaloulatlons  of  the  folloelng  type* 

1.  Lens  read  In 

2.  Third  and  fifth  order  oalculatlon 

3.  Automatic  correction  of  third  and  fifth  order  oalottlationa 

U*  Skew  ray  tracing  modes 

5.  Radial  energy  distribution 

6.  Frequency  response 

7«  Library  routine 

The  ORDEALS  Optics  Prof>raina 

A  aet  of  programs  have  been  prepared  to  carry  out  many  of  the 
requirements  for  automatic  lens  desl^  on  an  IBl  7070.  This  set  of 
programs  Is  oalled  ORDEII^. 

The  ORDE/UiS  programs  are  desombed  In  the  following  widteups  idiloh 
have  been  prepared  for  genered  circulation.  These  notes  were  prepared 
for  a  Rochester  Summer  School  where  the  programs  were  deeorlbed  and 
desunstrated.  The  programs  written  up  ai*ei 
1.  Optics  Programs  for  the  IBM  7070  Computer 

This  writeup  Is  available  in  limited  supply.  It  Is  also  being 
constantly  liiproved  and  added  to.  (Included  with  this  report  as  Appendix  A) 
n.  A  general  Llnearlytlon  Method  for  Automatic  Lens  Correction. 

This  paper  has  not  yet  been  presented  for  publloatlm  because  all  ■ 
the  latest  features  in  this  program  have  not  been  checked  out.  (Inoludad 
with  this  report  as  Appendix  B) 
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III»  Introdaotlon  to  the  OeCTnetrlc  Optical  Freqtiency  Response 

This  report  deeorlbes  the  theory  Justifying  the  use  of  Oeonetrloel 
fre<iuenoy  reeponee*  (Included  ulth  this  report  as  Appendix  C) 

Future  Work 

This  OROBAL  program  Is  now  being  put  into  subroutines  in  Fortran 
language  suitable  for  a  7090.  The  new  program  will  use  automatic 
oorreetlon  based  on  ray  deviations  rather  than  using  aberration  theory. 

Ey  tracing  13  rays^  it  is  possible  to  evaluate  the  optical  system  up 
to  the  7th  order.  After  experience  with  this  method  has  been  evaluated 
the  decision  will  be  made  as  to  whether  to  add  the  third  and  fifth  order 
aberration  theory.  The  latest  method  will  optimize  using  ge<»netrloal 
ray  aberration.  At  a  much  later  stage  it  will  become  necessary  to  in¬ 
corporate  wave  front  aberration  correction. 

Publications 

IXirlng  this  period  we  have  published  the  following  paper  in  the 
Journal  of  the  Optical  Society, 

Creative  Thinking  and  Computing  Machines  in  Optical  Design.  Robert  E. 
Hopkins  and  Gordon  Spencer.  J.  Opt.  Soc.  Am.  7ol.  No.  2,  172-176, 
February  1962.  This  paper  is  Included  in  this  report. 

The  ORDEALS  program  writeup  and  tape  units  have  been  made  available 
to  three  optical  companies. 

1.  Texas  Instruments 

2.  F.  M.  A. 

3.  Pacific  Optical  Company 
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1h«  Ibllevlng  oompanles  have  used  the  prograa  at  Rooheeter. 

1.  Paoiflo  Optical  Oo. 

2.  Batisoh  and  Lend: 

3.  Tropel  Ino. 

It.  Eastman  Kodak 

Itek  Corporation 

Personnel 

This  oontraot  has  supported  the  graduate  work  of  Cordon  Spencer. 

He  la  now  writing  a  PhD  thesis  on  this  problem  of  autosuttlo  correction. 

X 

The  oontraot  has  also  partially  supported  Mathew  Rlinner.  Mr.  RlaaBor 
Is  now  completing  a  Masters  thesis  on  Fifth  Order  Aberration  Theory. 

Mrs.  Alice  Sinclair  has  bewi  used  to  write  programs.  Her  main 
contribution  has  been  in  preparing  a  Library  routine.  This  feature  allows 
us  to  store  for  rea4y  aooesa  all  optical  design  solutions. 


TABU  Of  oowro 


nRMBOOlKMi 

A*  The  Ordeals  ProgrMS* 

1)  There  are  three  segMnts  of  prognM. 

B«  Basle  Featttrea  of  Call  and  Data  Oards. 

1)  first  word  alphabetlo. 

2)  ffpa  of  floating  point  igrataai. 

0.  MUB  u?  or  ni  uhs  m. 

1)  uMn. 

2)  OoBMant  cards. 

3)  Index  seleotlon  card. 

li)  Oorrature,  thickness  and  index  cards  for  eracy  sarfaoe 
including  the  object  and  ijsage. 

5)  The  paraxial  raj  data. 

6)  Call  cards  for  Other  ptiograu. 

0.  nmi  GALL  ADD  DATA  GAROB. 

1)  UBIIH. 

2)  Co— ant  oards. 

3)  Index  seleotlon  card. 

It)  Oarratore  oards. 

T^rpes  of  ourratore  oards. 

Prerlous  ourrature  pick  up. 

Conic  sections. 

Aspheric  surfaces. 


5)  SdataiMs  earda. 

of  thidoioaa  oarda* 

OOLoar  a^artnro  aolTO* 

Objoet  diatanoa. 
i)  Indi  of  rofvootlflB  oarda. 

a)  •afraotlng  aarfaoo. 

b)  lafloettng  aorfaeo. 

7)  Puwdal  laj  Data. 

S)  POU. 
f)  OHO. 

a)  fair  oarda. 

b)  OOni  oarda. 
a)  AID  oarda. 
d)  «>IJir  oarda. 
a)  mU)  oarda. 

f )  Itau  lAiioh  augr  bo  onrootad. 

g)  SfidiTldaal  aorfaeo  iXmm  idaieb 

10)  ADfD 

U)  lap  tboeo  Nodoa. 

a)  Axial  Fan. 

b)  foil  Fiald  Traoo. 

a)  XmUvldaal  Ragr  Traoo. 

1)  Tkwafor  oard  HDTR 

2)  Fraotional  ebjoot  hol^ 


It)  OlMor  ApMrior**. 

a)  AntoMtlo  Maignaant. 

b)  QLaar  apartfure  oarda. 

e)  OQTQA.  Oatalda  elaar  apnimft 

d)  IRGd.  Inaida  ^Uar  a|>annfa. 
13)  %•%  diagraa.  SAM 

a)  Wara  langth  «al#it. 

b)  firaetional  bbjaet  halgbt. 

e)  Fbona  Aift. 
d)  f^pot  Hat. 

lit)  ladlal  liargjr  Oiatrlbatlon. 

Oantar  dLaplaoaawnt. 

15)  Ffaqaanejr  Raaponaa. 

16)  Lana  Raad  Oat, 


B. 

OOTPUI. 

27 

P. 

SAMPLE  ntm  AMD  OUTPUT. 

31 

0. 

Altarlng  Fljcad  Oonatanta. 

31 

H. 

Oparatlon  of  tha  Maohlna, 

32 

I. 

Qanarol  CooaMnta, 

33 
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Appendix  A 

opncs  FsoaiuMS  idr  ihb  7070  ooNivm 
Migr  10,  1962. 

A  aet  of  design  and  sralttetlon  prograau  oellad  Ordeals  here  bean 
Witten  for  the  I.B.M.  7070  machine.  Ihe  following  writer  is  a  brief 
description  of  the  program  Input  and  output  features.  The  prograai  la 
being  constantly  worked  on,  so  this  description  la  subject  to  alteratBans. 
There  has  been  no  attempt  in  this  Hrlte*^p  to  provide  a  coaplate  dasori|ition 
of  the  program,  but  It  la  available  In  autocoder  language. 

There  are  236  pages  of  pn^ram.  So  far  there  are  7$  pages  of  flow 
diarts. 

The  ordeals  programs  are  written  for  an  I.B.M.  7070  machine  which 
eontains  the  following  basic  units. 

1)  Three,  channel  one,  tj^>e  units 

2)  Three,  channel  two,  tape  units 

3)  A  card  reader  unit 

U)  A  card  punch  unit 

$)  A  lltOl  printer  unit 

6)  A  lOK  storage  7070. 

The  programs  available  aret 

1)  Autosutlo  first,  third  and  fifth  order  aberrations  plw 
seventh  order  spherical  aberration, 

2)  A  ragr  trace  program, 

3)  Spot  diagram,  radial  energy  distribution  and  pseudo 


goosMtrioal  frequency  response 


The  pragram  oan  aooaimodAte  centered  rotatlonally  eymnetrlcal  eurfaoea 
of  the  fonn 

Z  -  ^ _ ♦  (e)3^  ♦  (f)s6  ♦  (g)S®  ♦  (h)S^O^  (1, 

1*  \/l  -  (K+1)C?s2 

ahere  3  -  ]/  ♦  l2. 

Oeoenterlng  and  tilt  pliia  ncai-rotatlonally  symetrioal  surfaces  are 
being  added  to  the  program.  A  write-up  of  Spxyte  (Special  Purpose  RAI  TRACI)  la 
included  In  Appendix  #1. 

The  7070  Ordeals  program  has  been  planned  and  %irltten  by  Mr.  Mathew 
Rlntier  and  Mr.  William  Hennessy.  Mr.  Rimmer  wrote  the  aberration  analyalSf 
and  Mr.  Hennessy  planned  and  wrote  most  of  the  other  features  In  the  program. 
Mr.  Gordon  Spencer  has  worked  out  the  details  of  most  of  the  computation 
techniques,  and  formulated  the  automatic  correction  method  used.  (See 
Appendix  2.)  Dr.  Robert  Hopkins  has  mostly  contributed  by  Insisting  that 
the  program  be  easy  to  use  and  require  a  minimum  of  input  data. 

A.  The  Ordeals  Programs. 

All  of  the  7070  Ordeals  programs  are  located  on  tape  10.  The  major 
programs  are  called  Into  core  storage  for  execution  with  call  cards.  The 
call  cards  are  simply  cards  with  the  name  of  the  program  punched  on  them. 

There  are  three  segments  of  programs.  One  set,  called  the  I  segment,  remains 
in  core  storage  once  it  has  been  called  in.  The  programs  in  segments  one 
and  two  have  to  share  space  in  core  storage.  If  one  calls  programs  from 
segment  two,  the  programs  from  segment  one  >hldi  are  in  core  are  wiped  out. 

The  programs  with  their  call  names  and  segments  are  listed  belowi 

Segment  I. 

LENIN  Lens  read  in. 

Lens  punch  out. 


LOAD 


Load  constants 
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Segment  1. 

AXTAX 

Axial  fan. 

rUFTR 

Full  field  traoe* 

miTR 

Individual  ray* 

8IM 

Spot  diagram. 

SHIFT 

Fbeus  shift. 

RED 

Radial  energy  distribution. 

SegRunt  2. 

FORD 

Fifth  order 

OORD 

Fifth  order  correct. 

AUTO 

Automatic  correction. 

Segnant  1  eontalna  programs  used  in  connection  with  ray  tracing.  Sag- 
Mnt  2  oontaina  programs  connected  to  aberration  calculation  and  correction* 
SepHnt  I  oontaina  programs  coiwaoo  to  ray  tracing  and  aberration  analysis* 

B*  BASIC  FEATURES  OF  CALL  AND  DATA  CARDS. 

Isput  cards  are  always  pundted  with  alphabetic  infoznatlon  in  the  first 
field.  Either  alphabetic  or  nuiaaric  Infonnatlon  may  be  punched  in  the  fire 
remaining  fields.  Fields  are  separated  on  the  card  by  one  or  more  blank 
eolumna.  The  first  field  may  be  punched  anywhere  on  the  card*  but  usually 
It  Is  started  In  the  first  column.  A  field  may  contain  a  of  tsn 

alphabetic  eharaotera  or  twenty  numeric  characters. 

Di  aoat  oases  y  numexrlcal  Information  is  pundied  In  floating  point  form, 
vltli  tho  sign  flivty  the  exponent  next  and  the  nuad>er  last.  Ibr  eocnplsi 
$11000000  represents  the  number  ♦  1.0  and 
-  Ii91000000  represents  the  noaber  •  0.01. 
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lh«  floating  point  exponent  refers  to  $0  m  a  baae,  and  exponants  majr  go 
from  99  to  01.  It  la  not  naoaaaarjr  to  punch  the  Inalgnlfioant  leroa  in  a 
nuaher}  henos  the  number  •O.OOUOOOOO  could  be  punched  and  the  above 

numbers  aa  $11  and  >lt91« 

If  all  the  paramstera  on  a  oard  are  aero,  the  alphabetic  field  la  the 
only  one  ihloh  must  be  pundied.  If  the  first  parameter  la  asrOf  but  the 
eeoond  la  not*  than  It  la  neoassasy  to  punch  a  single  aero  In  order  to  keep 
the  proper  order  of  the  wrds.  For  exaaqile,  an  aspherlo  oonstanta  oard  may 
have  four  fields  plus  the  alphabetic  field.  If  the  first,  second  nd  fourth 
cmiotants  are  aero  and  the  third  constant  la  -0.000026,  It  auqr  be  punched 
as  follows t 

ASffl  0  0  .U626. 

C.  MAKE  UP  OF  IHE  LENS  DECK. 

Since  there  are  several  options  available  for  setting  up  a  lens  deok, 
a  general  outline  la  presented  here. 

1.  LENIN.  This  call  card  calls  the  lens  read-ln  program,  reads 
the  lens  data  and  executes  a  paraxial  trace. 

2.  Ccranent  cards t  Cards  may  be  inseid;ed  to  identify  lane, 
give  date,  etc. 

3*  Index  selection  cardt  This  card  defines  the  wave  lengths 
to  be  used  in  the  third  and  fifth  order  aberration  cosqm- 
tatlons.  Five  wave  lengths  are  specified!  mid-wave  length, 
two  for  primary  color  aberrations  and  two  for  seooodaxy 
color  aberrations. 

U.  First  curvature  oardi  defines  the  curvature  of  the  object 


surface 
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$,  First  thloknass  osrdt  defines  the  axial  separation  betwean 
the  object  surface  and  the  first  surface  of  the  lens  system. 

6.  Index  of  Refraction  oardt  This  card  defines  five  indices 

of  refraction  (usually^  A,  C,  d,  F«  g)  for  the  medium  follow¬ 
ing  the  surface.  Any  five  wave  lengths  nay  be  used,  but 
they  are  always  referred  to  as  A>  0,  D',  F  and  Q  (In  that 
order)  in  the  program. 

7.  Ourvature,  thickness  and  Index  cards  for  the  rasudnlng 
surfaces.  Thickness  and  Index  cards  refer  to  the  space 
following  the  surface. 

8.  A  final  curvature  card  for  the  final  Image  surface. 

9.  The  paraxial  ray  data  card. 

10.  Call  and  data  caids  for  other  programs. 

D.  INPUT  DATA 

Details  of  Input  and  Call  Cards. 

1.  LEMINt  This  Is  a  call  card  for  'Uie  lens  read-ln  program. 
LENIN  may  be  punched  in  any  position  In  the  card.  All  the 
letters  must  be  together^  however,  with  no  spaces  between 
them. 

2.  Comment  cards  i  Following  the  LENIN  card,  one  may  add  as 
many  oomnent  cards  as  desired.  A  comment  card  must  stairt 
with  an  asterisk  In  woxd  field  one.  Comments  In  numerlo  or 
alphabetic  form  may  be  punched  according  to  the  following 
rules I 

a.  Words  should  not  exceed  ten  digits. 

b.  Do  not  extend  comments  beyond  column  59. 
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0.  There  must  be  a  space  between  the  asterlak  and  the 
first  diaraoter. 

A  card  with  an  asterisk  may  be  used  as  a  spacing  card* 

Sxaaqplet  The  coanents,  "John  Saiith,  Teleseope  ObJeotlTSf 
OeoaBd>er  8,  1962”,  could  be  pundisd  In  the  following  Hoar 
esrdat 

* 

«  JOHN  SMITH  TEUSCQPB  OBJECnVK 
*  DEC  8  1962 
« 

3«  Ihdex  of  Refraction  selection  cardt 

Field  one  la  punched  with  the  alpha  word  NRSEL.  Field  two  speci¬ 
fies  the  wave  lengths  used  In  the  calculations,  the  first  letter 
designating  the  major  color.  This  wave  length  will  be  used  for 
the  paraxial  trace,  third  and  fifth  order  aberrations.  The  next 
two  letters  refer  to  the  primary  chromatic  aberrations)  the  last 
tw9  letters  refer  to  the  secondary  chromatic  aberrations.  Any 
permutation  of  ACOFO  with  repetitions,  making  a  total  of  five 
letters.  Is  permitted.  Note  that '”FC”  and  "FD”  will  give  chromatic 
aberrations  with  the  usustl  slgn'sonventlons,  while  *0F”  and  ”DF” 
will  give  aberrations  mith treversed  signs. 

Example  t  To  cio  third  order  and  fifth  order  calculations  In  D 
light,  primary  color  between  F  and  C  and  secondary  color  between 
F  and  0,  punch  the  NRSQ,  card  as  follows  t 
NHSEL  SFGFD 

U*  Curvature  Cards i 


There  are  the  following  types  of  curvature  cards  which  nay  be 
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need  in  the  programt 

Card  Type 

Code  for  Field 

a)  Regular  curvature 

CF 

b)  Regular  curvature,  stop  surface 

CVS 

c)  Regular  curvature,  image  surface 

CVI 

d)  Angle  solve  on  the  axial  ray 

UA 

e)  Angle  solve  on  the  chief  ray 

UO 

f)  Radius  of  curvature 

RD 

g)  Previous  curvature  pick  up 

PC 

h)  Previous  surface  pick  up 

PS 

i)  Previous  radius  pick  up 

PR 

The  S  may  be  added  to  the  code  for  any  curvature  card  In  order 
to  make  the  surface  the  stop  surface.  For  exaaple,  PCS,  UA8  are 
allowed  codes  for  curvature  cards.  An  I  may  be  added  to  any  ounrao 
tore  code  in  order  to  make  the  surface  the  Image  surface,  e.g.,  PCI. 

The  first  curvature  card  in  a  lens  system  is  always  the  objeet 
ourvatui*e  card.  Usually  It  has  a  curvature  of  aero,  so  It  is 
pundied  CV.  The  zero  Is  not  needed. 

Previous  curvature  pick  up>  Often  it  is  necessary  or  desirable 
to  make  a  curvature  identical  with,  or  hare  a  constant  dlffer«Me 
from,  a  previous  curvature.  This  is  dme  with  the  code  PC.  In 
field  two  the  surface  number  of  the  surface  to  be  picked  up  is 
inserted  in  fixed  point.  The  cuirvature  may  be  picked  up  with  a 
reverse  in  sign  by  punching  a  minus  sign  in  front  of  the  surface 
number.  A  constant  difference  (t  C)  may  be  added  (or  subtracted) 
to  the  picked  up  curvature  by  inserting  a  floating  point  number 
in  field  three.  The  conic  constant  K  will  alwiqrs  be  picked  up,  and 
Its  sign  will  not  be  changed. 
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Ejcanplet  Siippose  that  va  vant  the  Ifth  surface  to  pick  up  the 
curvature  on  the  10th  surface  with  a  reverse  in  sign  and  a  curvature 
difference  of  •4.01.  Pundi  the  l^th  curvature  card  as  foUouet 

PC  -10  -li91. 

If  the  code  PR  ie  used,  it  picks  up  the  previous  radius  and 
it  adds  a  AR. 

If  the  code  is  PS,  it  picks  up  all  the  surface  constants  of 
the  previous  surface,  including  the  constant  K  (see  Equation  1), 
and  the  asiherio  constants  if  there  are  any. 

Conic  sections t  Any  floating  point  numbers  inserted  in  the  third 
word  field  of  a  CV  card  are  Interpreted  as  the  conic  section 
constant  K.  The  following  table  relates  the  conic  section  K  to 
the  various  types  of  conic  sections  available t 


Type  of  Conic  Section  Constant  K 

Paraboloid  K  ■  -1 

Hyperboloid  K  ^  -1 

Ellipsoid  revolved  about  the  major  axis  -1<K<0 

Ellipsoid  revolved  about  the  minor  axis  K  0 


p 

The  conic  section  K  ••  -E  where  E  is  the  eccentricity. 

Exanplet  A  paraboloid  could  be  described  as  CV  -I49I  -511. 

Aspheric  Surfaces t  General  aspheric  surfaces  of  the  type  given  in 
aquation  1  may  be  Inserted  at  any  time  by  adding  a  card  directly 
after  a  curvature  card^  with  ASFH  in  field  one.  Fields  two,  three, 
four  and  five  are  reserved  for  the  coefficients  e,  f,  g  and  h. 

(See  equation  1.) 
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An  aspheric  surface  idth  defonMtlon  from  a  paraboloid 
would  than  have  two  cards  as  follower 

CV  -U91  -511 
ASFH  0  0  -1)626  •‘Uiih 

5*  Ihlckness  cards i 

There  are  several  types  of  thickness  cards,  as  Indloatad  in 
'  the  ftollowing  tablet 

Card  Type  Ooda  for  Flald  1 

a)  Regular  thickness  1R 

b)  TA  solve.  The  computer  conpites  the  Ik 

required  thickness  to  make  the  y  of 

the  axial  ray  on  the  next  surface 
have  the  prescribed  YA. 

o)  YG  is  the  sme  for  the  dtief  ray.  TO 

d)  Previous  thickness  pick  up.  This  FT 

is  exactly  analogous  to  the  previous 
curvature  code. 

e)  Clear  aperture  solve.  CA 

Olear  apertuM  solver  This  feature  provides  for  appropriate  axial 
separations  between  surfaces.  There  are  two  types  of  adjaoaiit 
surfaces.  They  are  called  closed  if  they  edge  contact  when  brou^t 
together.  If  they  contact  at  the  center,  they  are  called  open. 

If  the  code  CA  is  used  cm  the  j  surface,  the  progran  oheoks  to 
see  what  the  (j  1)  surface  is.  Next  it  oheoks  to  see  whether 
the  material  between  the  two  surfaces  is  glass  or  air.  Finally 
it  oo^putes  the  thickness  between  the  ^  and  the  e  1}  surfaoes 
aooordlng  to  the  following  sdiedulet 
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Oondltlon 

Air 

Olaas 

)(  Open 

THj-a 

a  noeilnally  la  0.01 

IHj  -  b  CA. 
b  nominally  Is  0.2 

0  Closed 

-*3  .  1 

™  -  *J  -*J  ♦  1  *  c 

0  nominally  Is  0.1 

In  no  OMo  la  tha  glaaa  thloknaaa  nada  loaa  than  0.2  CA.  lha  aooond 
flald  of  tha  GA  card  oontalna  tha  nunarloal  Talue  of  tha  radlua  of 
tha  olaar  apai^ura  raqtilrad. 

Fbr  axaafilat  A  clear  apartura  radlua  of  1.2  would  be  punohadt 

CA  5112. 

lha  aaga  Zj  and  Zj  ^  1  computed  with  the  following  approxl- 
nata  formulat 

2-  032  ^1  zC2jj2L(K  +  D)  /  (2  -C2  S2  (gel)  . 

Qha  ahould  remember  that  the  edge  thlokneeaea  (c)  of  the  lenaas 
are  datarmlned  by-  tha  fixed  conatant  0.1.  Thla  conatant  lugr  be 
altered  at  any  tine,  aa  well  aa  the  conatanta  a  and  b,  by  Inaerting 
LOAD  cards  after  tha  fORD  card.  See  the  section  (O)  on  altering 
program  constants. 

Object  ■  Pjatanoet  In  this  program  all  object  distances  THg  are 
considered  to  be  finite.  To  make  an  infinite  object  distance, 
anm  merely  needs  to  use  a  large  number  for  TH,^.  For  example.  It 
is  otmranlant  to  use  a  floating  point  number  such  as  601,  or  60F, 
where  F  is  tha  focal  length  of  the  system.  This  is  useful  because 
than  tha  object  halgiit  has  the  same  value  as  the  image  helf^t, 
with  merely  a  change  of  the  exponent.  For  example,  suppose  lens 
has  a  focal  length  of  2U"  and  It  is  to  cover  a  half  image  field  of 
6”.  Then  If  TH^  Is  made  602U,  Iq  Is  596.  There  Is  one  danger  with 
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this  If  th«  first  surfaoa  la  piano  or  allgbtljr  ooneatvo  towards  tha 
objaot.  Than  which  Is  aqaal  to  yo  ♦  u.^  aajr  baoena  so  snail 
that,  whsn  It  Is  ralsad  to  tha  7th  powsr,  as  Is  dona  in  ealoalatim 
tha  7th  ordar  spharloal  abarratlon,  It  asgr  undarflow*  If  this 
happans,  ona  may  raduoa 
6.  Index  of  Rafractlon  Cards t 

a)  Bafraotlng  surfaoet  Punch  INDEX  In  word  flald  ona  and  flra 
floating  point  Indices  in  word  fields  two,  three,  four,  flra  and 
slat. 

Tha  five  Indices  are  referred  to  liqr  tha  alphabetic  letters 
A,  Qy  0,  F  and  0.  One  may  insex^  any  indices  desired,  but  they 
Bust  be  referred  to  with  these  letters. 

Caution It 

1)  Punch  an  Index  other  than  zero  for  all  five  fields  evan 
though  they  may  not  be  used. 

2)  Ihe  index  card  Is  a  potential  source  of  grave  danger. 

A  mlspunched  number  hex^  can  mean  a  false  index  of 
refraction  and  go  unnoticed  until  the  final  design  Is 
made  up.  Be  sure  and  check  these  Indices!  It  Is  advls- 
abln  to  compute  V  numbers  directly  from  the  numbers  on 
the  card.  It  Is  also  advisable  to  make  up  a  set  of 
prepunched  cards  of  the  comnonly-oised  glasses. 

b)  Reflecting  surfacest  If  a  surface  la  to  be  reflecting,  simply 
punch  REFL  in  word  field  one  In  place  of  INDEX.  No  other  nundtars 
are  needed.  The  program  automatically  Inserts  negative  index 
values  for  all  the  following  surfaces  until  another  REFL  card  Is 
encountered.  Daring  this  rlgh^-to-left  travel,  one  must  Insert 
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thloknesaes  with  negative  values.  This  Is  me  place  where  the 
negative  previous  thickness  Is  useful. 

7*  Paraxial  Ray  Data  Cardt 

A  card  with  PXRAY  in  field  one  Is  used  to  provide  the  specl- 
flcatlon  of  the  axial  and  the  chief  ray. 

The  ItsBUi  entered  In  fields  two,  three  and  four  aret 
-  the  halt^t  of  the  object  at  the  object  surface.  It 
la  usually  chosen  to  be  e<iual  to  the  value 

to  be  used  by  the  lens. 

y^^  -  the  height  of  the  axial  ray  at  the  first  surface  of 
the  system.  It  Is  usually  chosen  to  be  the  maxi, mum 
height  acceptable  by  the  lens  for  axial  rays. 

y^^  -  the  height  of  the  chief  ray  on  the  first  surface. 

If  a  surface  in  the  system  Is  labelled  as  the  stop  surface, 

It  la  not  necessary  to  insert  a  value  of  y^.  The  program  will 
assume  a  value  of  3r^  and  then  trace  the  chief  ray  through  to  the 
stop  surface.  If  the  y  is  not  zero  on  the  stop  surface,  the  program 
autoiaatlcally  changes  y^  to  make  It  zero.  If  no  surface  Is  labelled 
stop  surface  and  y^  Is  left  zero,  then  the  first  surface  becomes 
the  atop  surface. 

8.  FORDi 

FORD  is  a  call  card  idilch  calls  in  the  aberration  correcting 
program  and  Initiates  the  computation.  Upon  coiiq>letlon,  the  progrm 
automatically  transfers  to  (X^. 

9.  OOROi- 

This  program  Initiates  the  automatic  correcting  part  of  the 
program.  The  first  thing  It  does  Is  ask  to  read  VARI  cards,  and 
finally  OORR,  ADD,  HOLD  AND  HOU)T  cards.  Those  cards  msy  be  oxplalnod 
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u  followit 

a)  VARI  cards. 

VABX,  cards  are  pundied  elth  ?ABX  in  field  one,  a  fixed  point 
Bttrfaoe  nunber  in  field  teo,  the  curvature,  oonlo  constant  k, 
aaphorlo  constants  -  or  thickness  to  be  varied  in  field  three,  wd 
a  fixed  point  veight  in  field  four.  Ihe  weights  vary  fltna  aero 
to  nine. 

fijcaml^t  To  vary  the  curvature  on  the  6th  surface  with  a  weli^t 
of  two,  punch  a  card  as  follows! 

VARI  6  CV  2 

It  is  possible  to  vazy  C7,  GC*,  TH,  AD,  AE,  AF,  AO,  BD.  If  a  YABX 
CV  card  is  used,  the  curvature  on  the  surface  will  be  varied  with¬ 
out  regard  to  the  curvature  code  used.  For  exanple,  if  the  surface 
has  a  UA  code  and  a  VARY  CV  card  is  used,  the  UA  will  be  varied. 

b)  CORR  cards. 

These  cards  specify  the  items  for  correction.  OORR  is  punched 
in  field  one.  Field  two  is  the  aberration  to  be  corrected.  Field 
three  la  the  target  value  for  the  correction.  The  itssM  idiiah  awgr 
be  corrected  are  listed  in  part  f  of  this  section. 

Isawplei  If  one  wishes  to  correct  the  spherical  aberration  8A5 
to  aero,  the  card  can  be  pun(died 

Qtm  3A$ 

o)  ADD  cards. 

If  a  (XRUl  card  is  changed  to  an  ADD,  than  the  aberration 
signified  in  word  field  two  is  added  to  the  previous  OORR  eard. 
never  pundi  target  values  or  weights  on  an  ADD  card. 


»  00  is  the  code  for  the  oonlo  section  constant  K 
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teanplet  One  may  correct  3A3  *  SA^  *  Ski  to  lero  by  Ineex'tlng 
the  following  three  cards t 

CORR  SA3 
ADD  SA5 
ADD  SA7 

One  could  correct  a  system  to  an  over-all  length  of  10*0  by 
inserting  the  following  set  of  cards. 

CORR  IH  1  521 

ADD  TH  2 

ADD  TO  3 

etc.  to 
ADD  TO  K-1 

d)  HOLDT  card. 

If  the  CORR  is  replaced  by  HOLDT  in  word  field  one,  the 
program  tries  to  minimize  the  difference  between  the  target  value 
and  present  value.  Field  two  contains  the  aberration,  field 
three  the  target  value  and  field  four  a  fixed  point  weight  which 
varies  from  zero  to  nine. 

Example  I  If  one  wants  to  minimize  the  difference  between  SA3  * 

SA5  *  Ski  and  aero  with  a  weight  of  2,  one  can  punch  the  following 
throe  cards t 

HOLDT  SA3  0  2 
ADD  SA5 
ADO  SA7 

The  HOLDT  cards  are  used  for  aberrations  which  may  not  land  thom- 
■ol^o*  for  complete  correction  to  a  target  value.  One  also  nay 

use  HOLDT  cards  to  correct  aberrations  after  all  the  variables 
are  used  up.  (hie  may  Insert  more  HOLDT  cards  than  variables. 
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•)  BOLD  earda. 

If  BOLDT  la  raplaoad  hf  BOLD  In  flald  eoa,  It  ImUeatM  Hint 
tha  praaant  ralua  of  tha  aberration  riioiild  baoona  tho  targot  anlM* 
Othorwlaa  tha  oard  la  axaotly  tha  aana  aa  tha  BOLD?  oard* 
f )  lha  foUowliig  ItaaM  nagr  bo  oorraotod  with  OOBB,  BOUT*  UD 


or  BOIi)  cardat 

rOGAL 

The  focal  length  of  tho  ayatcn. 

TAGH 

Primary  axial  color. 

TCH 

Primary  lateral  color. 

TACH2 

Secondary  axial  color. 

TCH2 

Secondary  lateral  color. 

SA3 

-  Third  order  apherlcal  aberration.  * 

aA5 

Fifth  order  aphezdcal  aberration. 

SA? 

Seventh  order  spherical  aberration. 

OOHA3 

Third  order  tangential  ocawi. 

0CMA5 

Fifth  order  linear  ctana. 

LOGHA 

Seventh  order  cubic  coma. 

TAS3 

Third  order  tangential  astigmatic  blur. 

TAS5 

Fifth  order  tangential  astigmatic  blur. 

DIST3 

Thix*d  order  distortion. 

0IST5 

Fifth  order  distortion. 

lOBSA 

Tangential  oblique  spherical  aberration 

aOBSA 

Sagittal  oblique  spherical  aberration. 

PTZ3 

Third  order  Petival  blur. 

PTZ5 

Fifth  order  Peti  blur. 

aAS3 

Third  order  Sagittal  astigmatic  blur. 

aA85 

Fifth  order  Sagittal  astipwtle  blur. 
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If  onlj  third'order  aberrations  are  asked  for,  the  pxogrw  does 
not  oosqpfttte  all  the  fifth  order  aberrations  during  the  antosuitlo 
oorreetlng  phase t 

g)  It  Is  possible  to  oorreot  sane  Itesu  on  IndlrLdual  surfaoes. 

Xn  this  oase,  the  OORR,  HOU),  MW,  or  H0IJ3T  are  punshed  as  fellSMSi 
Held  1.  CORR,  ADD,  HOLD  or  HOIDT. 

Field  2.  CODE  OF  SURFACE  ITBl. 

Field  3.  Surface  number. 

Field  I4.  Target  value. 

The  oodes  for  the  individual  surface  itesis  arei 
PI  Angle  of  Incidence  of  axial  ray. 

•  PIC  Angle  of  incidence  of  chief  ngr* 

PO  Angle  u  of  axial  ray. 

FUC  Angle  u  of  chief  ray. 

TH  Thickness. 

CORD  will  transfer  to  read  in  correction  data.  After  CORD,  one 
should  have  VARI,  OORR,  ADO,  HOIi),  or  HOLDT  cards.  These  oards  are 
followed  by  an  AUTO  card  described  in  the  next  section.  FORD  auto¬ 
matically  transfers  to  CUU)  upon  oonpletlon  of  the  fifth  order  prograsi. 
Thus  It  Is  not  neoessaxy  to  follow  a  FORD  card  with  a  CORD  card.  If 
FORD  Is  followed  by  a  call  card  for  another  program  Instead  of  a  TART, 
OORR,  ADD,  HOLDT  or  HOLD,  the  new  program  Is  called  and  control  la 
transferred  to  it.  CORD  may  be  \i8ed  to  oorreot  a  Ians  in  steps.  For 
mxjmplt,  after  FORD  one  mi^t  Insert  VARI  oards,  and  OORR  oards  to 
oorreot  axial  and  lateral  color  and  Petsval.  These  most  be  followed 
by  an  AUTO  card.  Then  a  CORD  card  will  permit  trading  in  awnobhro- 

Mtio  correction  data  and  an  AUTO  oard,  which  will  begin  oorreotloh 
froai  the  already  oolor-oorreoted  system. 
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lA*  ADTO  oalia  in  tha  aatoautio  oorraotlon  prograa  and  tmafara  to 
ialtlata  tha  autenatlo  eorraotlon  prooedora.  It  anat  ba  aaad  to 
laitlata  thla  faatora. 


Iha  autoMtlo  oorraotlon  nathod  la  briafljr  tha  foUowiiigi 
Tha  0,  t  tabla  idll  sTatanatloally  ba  TaTlad  ona  paraaatar 
at  a  tiaa.  Oapandlng  on  idiat  tha  daalgnar  aaka  the  aaohlaa  to 


oorraot,  it  vlll  ooafiata  third  and  fifth  order  abarratlana*  Will 
thla  rotttlna  It  la  poaalbla  to  ooaqmta  finite  dlffoMnoa  ratioa* 
After  all  tha  raqolred  dlffaranoa  ratloa  are  nna^mtad. 
are  loaded  into  the  aquation  aolTing  routine. 

Tha  aquation  solving  routine  aaaka  to  alnlalBa  a  fnnetlon 


aubjaot  to  the  alnultanaous  solution  of  the  set  of  aqnatlona» 

N 

In  Matrix  notation  this  leads  to  the  solution  of  the  folloidng 
aat  of  N  ^  L  equations  with  N  e  L  unknowns! 

a'^Ax  +BT\  «  A'^d  , 

Bx  -t-OX  -  e 

'^Mn/  \1>L| 
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Ihese  equations  are  used  in  the  design  program  to  solve  the  sets 
of  non-linear  equations  encountered  in  lens  design  bgr  assuming 
linearity  over  a  small  region. 

The  elements  of  the  x  matrix  become  pj  -  dlffarmaoe 

ratios  are  evaluated  at  ^^pj.  pj  is  the  value  of  Pj  at  a  predicted 
solution. 

The  elements  of  the  d  matrix  become  d^,  a  set  of  dlstanoes 
from  target  values  for  the  set  of  functions  to  be  minimised. 

The  elements  of  the  e  matrix  become  e^  \dilch  are  the  dlstanoes 
from  target  values  for  the  functions  to  be  corrected  exactly. 

Die  A  matrix,  is  a  set  of  Lagrange  multipliers. 

Ihls  method  is  the  general  form  for  McCarthy's,  Rosen's  mad 
MSimne's  method. 

A  more  detailed  analysis  of  this  procedure  is  described  by 
Cordon  Spencer  in  a  paper  presented  to  the  Optical  Society  of 
America  and  Included  in  Appendix  2. 

11.  Ray  Trace  Modes. 

There  are  four  modes  of  ray  tracing  available.  The  first 
three  modes  are  called  with  the  cards  AXFAN,  FUFTR,  INDVR,  The 
fourth  mode  is  SAN,  a  ray  tracing  for  spot  diagrams.  This  will 
be  described  in  Section  13.  The  first  two  call  in  programs  idiloh 
use  the  paraxial  ray  data  and  automatically  trace  through  a  group 
of  rays.  The  last  cairi  calls  for  data  specifying  individual  ray 
data.  A  descz*lption  of  these  modes  is  as  follows  i 

a)  Axial  Fani  This  mode  of  trace  is  called  by  a  card  with 
the  first  field  punched  AXFAN.  Immediately  following  this  card, 
one  must  have  a  card  punched  with  RADAT  in  field  one  and  the  wave 
length  code  for  up  to  three  wave  lengths  in  field  two. 
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BttHplet  To  trace  axial  fane  in  D  llg)it  and  0  lll^t*  VM  the 
foUoNing  two  oardat 

AXFAN 
RADAT  DC 


1b  do  More  than  three  wave  lengthif  extra  AXFAX  and  RABAT  barde 
are  needed. 

Ihe  program  traoea  three  ragre  In  the  flrat  wave  langth  from 
the  axial  object  point  (Tg  -  0)  at  the  fblloidng  valnea  of  X  on 
the  antranee  pupil  (1^^  ■  0). 


Hay  1  Xg„  -  0.5 

Ray  2  1^  •  0.75  yen 

Ray  3  Xen  -  1*0  y^ 

It  then  repeata  the  trace  for  the  other  wave  langtha  on  the  RADAT 
card. 

b)  Full  Field  Trace i  Thla  mode  of  ray  trace  la  called  with  a  card 
pundied  with  FUFTR  In  field  one.  It  la  immediately  followed  with 
a  RADAT  card.  Thla  RADAT  card  la  the  amae  aa  the  one  uaed  after 
AZFAMi  but  there  are  extra  fielda.  In  field  threo«  a  fixed  point 
number  apecifiea  the  number  of  off-axia  image  pointa.  Ihe  fixed 
point  nondMr  In  field  four  apeolflea  the  number  of  meridional 
raya  to  be  traced  on  each  aide  of  the  chief  ray. 
ftweplei  The  following  two  carda  are  a  oommon  oxMg>lo  of  the 


miRt 


FUFTR 

RADAT  DFO  2  3 
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Ihase  oarda  call  Iha  following  aequenoe  of  ragr  traolngi 
Ragra  1-3  An  axial  fan  in  D  ligfit 
”  An  aidal  fan  in  F  ll^t 

"  6-9  An  axial  fan  in  C  lic^t 

Ragr  10  A  ohiaf  ragr  in  D  li^t  fro*  an  object 

point  at  2/3 

*  11  A  chief  ragr  in  F  li^t  from  an  object 

point  at  2/3  y^. 

*  12  A  chief  ray  in  C  lifijit  from  an  object 

point  at  2/3  yo* 

Ragra  13**16  Rays  from  an  object  at  2/3  y^  but  confined 
to  the  meridional  plane  X  0.  The  ragra 
are  evenly  spaced  above  and  below  the 
chief  ray.  The  two  rays  fajrthest  frcai 
the  chief  ray  pass  through  the  top  and 
bottom  of  the  entrance  pupil. 

■  19-21  Three  skew  rays  from  the  object  at  2/3  y^ 

idilch  pass  thirough  the  entrance  pupil  with 
•  0.  The  rays  are  evenly  spaced  fron 

the  center  of  the  pupil  to  the  edge. 

"  22-2U  Chief  rays  In  0,  F  and  C  light  Itxm  an 

object  point  at  y^. 

*  25-30  Meridional  rays  like  rays  13-18  from  an 

object  point  at  y^. 

*  31-33  Three  skew  rays  like  ragra  19-21  from  an 

object  point  at  y^. 


Appendix  A 


-22- 


o)  IndlYlda*!  Ray  Tracei  For  Individual  nQr  trace,  one  auat  have 
the  following  arrangement  of  cardet 

1)  The  first  card  Is  a  transfer  card  with  IMDTR  punched  in 
field  one. 

2)  The  next  card  must  be  a  fractional  height  card.  FOBJH 
is  Inserted  in  field  one}  next  the  fractional  object  hel^t 
in  floating  point  is  placed  in  field  two.  The  rest  of  the 
card  magr  be  blank,  in  viiich  case  a  minimum  of  printing  will 
be  done  (no  print),  and  color  D  will  be  selected  for  the 
<diief  ray.  If  a  color  other  than  D  is  desired  (A  C  F  0) 

Oie  of  these  letters  must  be  punched.  If  desired,  PRINT 
may  be  punched  first  and  then  the  cblor. 

If  the  card  read  at  this  time  is  not  a  FOBJH  card,  the 
card  will  be  ignored  and  reading  will  continue  until  a  FOBJH 
c»rd  is  read,  or  a  transfer  card  to  the  next  routine. 

3)  The  next  cazd  must  be  EHPOO,  an  entrance  pupil  coordinate 
card.  X  on  the  entrance  pupil  is  punched  in  field  two  and 

T  on  the  entrance  pupil  is  punched  in  field  three.  The  rest 
of  this  card  may  be  blank  or  contain  a  color  (A  C  D  F  0) 
or  the  word  PRINT  or  the  word  NOPRT.  If  the  card  is  blank, 
printing  and  color  will  bo  controlled  by  the  FOBJH  card,  or 
if  that  was  blank,  no  print  and  D  will  bo  used.  The  print 
and  color  codes  are  originally  set  to  no  and  D  respectively. 
These  codes  can  be  changed  on  either  FOBJH  or  £NP(X)  cards 
and  remain  changed  until  a  new  change  is  read  in.  The  color 
and  print  words  can  be  entered  in  either  order,  or  either 
may  be  blank,  i.e,,  not  put  in  at  all. 
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12.  Clear  ApertvireBt 

Provlelon  is  made  in  'the  program  'to  stop  rays  ahldi  exoeed 
a  certain  clear  aperture  radius.  This  Is  done  to  out  dowi  on  the 
amount  of  useless  raytracing.  The  clear  aperture  radlua  may  be 
arrived  at  In  three  different  ways. 

a)  Automatic  Assignment!  If  the  designer  makes  no  attempt  to 
assign  a  clear  aperture  radius  to  a  surface,  the  clear  aperture 
la  autcmatloally  conputed  from  paraxial  ray  data.  The  clear 
aperture  Is  computed  from  the  following  formulai 

“3 -hi 

In  the  early  stages  of  design,  the  clear  apertures  are  assigned 
in  this  manner. 

b)  By  inseirting  a  CA  card  in  place  of  a  TH  cartl. 

c)  By  Inserting  an  OUTCA  card.  OUTCA  cards  are  Inserted  just 
before  the  thickness  cards.  To  provide  for  an  outside  clear 
^>erture  of  1.2.  the  card  would  be  punched  OUTCA  5112. 

An  OUTCA  card  should  not  be  used  on  a  surface  with  CA. 

The  OUTCA  card  Is  used  ^en  the  designer  wishes  'to  alter  the 
clear  cpeirturcs  assigned  automatically.  There  are  cases  where 
the  paraxial  ray  data  used  in  the  au'tomatlc  method  Is  not  accurate 
enough  and  rays  are  needlessly  blocked.  One  example  of  this  Is 
the  use  of  a  field  flattener.  The  formula  assigns  a  value  close 
to  CA  ■  y  because  y  is  very  small.  If  there  is  any  positive 
dlstox^lon.  alT  the  actual  rays  may  exceed  this  value  of  y.  For 
'this  case  It  is  advisable  to  insert  an  OUTCA  card  on  the  field 
flattener  surfaces  in  order  to  overrule  the  automatic  aaslgnment 
of  CA.  0UTC4  cards  are  used  extensively  in  the  later  stages  of 
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•  design  idien  the  energy  distribution  is  computed* 
d)  Die  pngram  also  provides  for  a  card  called  INCA.  Dlls  card 
blocks  out  all  rays  closer  to  the  optical  axis  than  a  specified 
Inside  clear  aperture  radius.  If  a  surface  has  both  an  ODTQA 
and  INGA,  they  are  inserted  on  two  cards  Immediately  preceding 
the  thickness  card.  Dielr  order  does  not  matter. 

Die  INGA  card  does  not  block  out  chief  rays,  since  thoM 
rays  are  used  as  references. 

13.  Spot  Diagram t 

SAM  Is  the  call  card  for  the  spot  diagram  program. 

a)  SAM  should  be  followed  immediately  by  a  card  with  WAVEM  In 
field  one.  Die  fields  two  to  six  In  this  card  are  reserved  for 
fixed  point  weights  for  the  five  wave  lengths  on  the  Index  of 
refraction  cards.  Die  weights  are  used  to  specify  the  nundier 

of  rays  to  be  traced.  Die  weight  times  ten  Is  the  number  of  rays 
traced.  For  example!  a  UAVEW  card  punched  1  0  3  0  1  would 
trace  10  rays  In  A  light,  30  rays  in  D  light  and  10  rays  In  0 
light.  Die  spot  diagram  contains  coordinates  for  all  the  waive 
lengths . 

b)  Immediately  following  WAVEW  Is  a  card  with  FOBJH  In  field  one. 
Field  two  Is  resez^red  for  the  floating  point  fractional  object 
height.  Dlls  card,  therefore.  Indicates  the  object  height  for 
the  spot  diagram.  It  Is  necessary  to  Inseirt  a  FOBJH  card  for 
each  object  height. 

o)  If  a  focus  shift  Is  needed,  a  card  with  SHIFT  Is  Inserted  after 
the  FOBJH  card.  Die  second  field  on  the  SHIFT  card  should  contain 
the  focus  shift  as  a  floating  point  number. 
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d)  The  spot  diagram  ooordlnatas  may  be  listed  by  Inserting  a 
card  with  SPLIT  in  field  one.  All  the  preceding  spot  dlagrms 
will  then  be  listed. 

The  spot  diagram  program  contains  a  subroutine  that  uses 
the  paraxial  ray  data  and  the  clear  apei*tures  of  the  lenses  and 
computes  the  shipe  and  area  of  the  vignetted  aperture  on  the 
entrance  pupil.  The  rays  that  are  traced  for  each  color  are 
evenly,  distributed  over  this  area, 
lii.  Radial  Energy  Diatributioni 

A  call  card  labelled  RED  calls  in  the  program  to  compute 
the  radial  or  encircled  energy  distribution  for  all  the  preceding 
fractional  object  heights.  This  program  computes  the  per  cent 
radial  energy  distribution  at  five  per  cent  intervals.  The  circles 
are  all  assumed  to  be  centered  .iround  the  chief  ray  in  the  major 
color  as  specified  on  the  NRSEL  card.  It  is  possible  to  shift  the 
center  of  the  circles  by  Inserting  a  card  with  CD  in  field  one, 
followed  by  the  floating  point  displacement  in  fields  two  and  three 
for  X  and  Y  center  displacement  respectively.  CD  stands  for  center 
displacement.  One  may  insert  as  many  center  displacement  cards 
as  needed. 

15.  FRBQUBICY  RESPONSE 

This  program  computes  a  geometidcal  frequency  response  from 
energy  distribution  data.  The  formula  used  is 
T  (w)  -  J.Ccjr) 

f 

wiiere  E  is  the  difference  between  two  successive  values  of  the 
energy,  E.  r  •>  (R)(RSC),  where  R  is  the  arithmetic  mean  of  the 
two  E  values,  and  RSC  ■  (sec  idiere  0  is  the  obliquity  angle. 
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W  •  2tfH  where  N  ie  the  number  of  epatlal  oyclee  per  unit  length* 
R  end  N  must  be  euoh  that  wr  la  dlmenelonleeB*  Division  by 
nonnallsee  the  response  function* 

There  are  two  entries  to  the  programs 


1)  FREQT* 

This  call  card  will  initialize  coiiq>utation  of  energy 
distributions  for  each  spot  diagram  computed  thus  far  on 
the  present  lens  sdiich  will  appear  in  the  output* 

a)  The  obliquity  factor  and  up  to  twenty  values  of  N 


are  read  in  as  follows i 

RSC  *$11 
FRQ  +511 
PRQ  +512 
FRQ  513 
FRQ  531 


b)  The  output  is  two  columns, 
and  the  other  the  response 


one  showing  the  frequanoy, 
at  that  frequency* 


2)  FREQC* 

This  card  calls  for  energy  dj.stribution  data  to  be 
read  from  up  to  twenty  cards  with  the  code  DIST  in  field 
one,  the  energy  in  field  two  and  the  radius  in  field  threes 


DIST 

1495 

I485 

DIST 

501 

I489 

DIST 

511 

1492 

IMS 
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Die  progrsm  then  transfers  to  read  In  the  obliquity 
factor  and  the  frequencies,  as  In  lA. 

16.  Lens  Read  Outi 

A  call  card  marked  LEO  Is  used  to  provide  an  updated  lena 
deck.  It  may  be  inserted  any  place  In  a  deck,  provldod  It  doean't 
separate  a  call  card  and  Its  required  data  cards. 

E.  OUTPUT. 

The  output  format  for  the  7070  programs  are  ahoun  on  pages  29-30. 

The  first  column  Is  the  line  list  number*  The  second  oolunn  of  three 
numbers  indicates  the  type  of  printing.  The  ^00  line  is  the  ln|mt  data 
obtained  from  the  PXRAY  card  and  the  object  distance  THg* 

The  lines  $01  to  506  are  the  lines  printed  for  each  surface  in  the 
lens.  The  second  column,  which  is  a  series  of  01  values,  is  the  surface 
number.  In  this  example  only  the  first  surface  is  shown.  The  aberra¬ 
tions  shown  are  the  surface  contributions. 

The  lines  507  to  510  are  the  total  third  and  fifth  order  aberrations 
for  the  system. 

The  line  511  shows  the  printing  during  each  change  cycle.  Vq  Is 
the  initial  value  of  a  parameter.  is  the  new  value.  Vi  -  Vq  Is  the 
change  in  the  parameter. 

The  lines  200  to  253  are  the  lines  indicating  the  printing  in  ray 
tracing.  Hie  last  digit  in  the  number  indicates  the  wave  length, 
according  to  the  following  rulet 

A  -  1 
C  -  2 

D  -  3 
r  -  h 
0-5 


Lin*  223  shows  the  full  print  on  individual  surfaces. 

the  first  line  of  printing  after  RED  has  the  1^0  printing  oode. 

The  line  of  printing  Is  shown  on  page  30.  The  Itesis  printed  In 
this  line  arei 

1.  nie  area  of  the  vignetted  aperture. 

2.  R^  la  the  radius  of  the  aperture  Ujmltlng  the  upper 
rays,  as  projected  on  the  entrance  pupil  plane. 

3.  Is  the  center  of  this  aperture  as  located  In  the 
entrance  pupil. 

U.  R]^  Is  the  radius  of  the  aperture  limiting  the  lower 
rsya,  as  projected  on  the  entrance  pupil  plane. 

5*  la  the  center  of  the  lower  aperture  aa  located  in 
the  entrance  pupil. 

6.  Rq,)  Is  the  radius  of  the  entrance  pupil. 

Following  this  1^0  line,  the  per  cent  energy  table  Is  printed. 

Above  the  table  four  numbers  are  printed,  and  labelled.  Ihey  aret 

1.  FOBJH  fractional  object  height. 

2.  Focus  shift. 

3.  CD|,(  center  displacement  In  x  direction. 

1(.  CD|y  center  displacement  In  y  direction. 

Ihe  per  cent  energy  table  Is  made  up  of  two  columns,  the  per  cent  energy 
and  the  corresponding  radius  of  spot. 

Following  the  table,  a  single  line  of  fixed  point  numbers  gives 
the  number  of  rays  traced  In  each  color  and  the  total  number  of  rays. 


7070  ODTPDT 


fH 

O 

o 

O 

iH 

O 

rH 

O 

iH 

O 

VTV 

CM 

O 

u\ 

O 

« 

vr\ 

507 

8 

\A 

& 

U\ 

V\ 

CHANGES  IN  V. 
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F.  SAMPLE  INPUT  AND  OUTPUT. 

A  sanpla  Output  from  FORD  Is  shown  on  pages  2  and  3  In  Sonpla  Design 
number  three.  Tlie  program  first  pj-ints  out  the  Input  data  from  LBflN  to 
PXRAY. 

On  Page  ij  the  Input  from  FORD  to  the  final  CORR  DIST3  is  shown.  The 
remaining  lines  from  0U2  to  05^  show  the  changes  made  In  the  first  iteration 
and  the  total  aberrations  for  the  first  iteration. 

Page  15  shows  the  input  cards  AUTO,  LEO  and  FUFTR.  Tlio  remainder  of 
page  15,  pages  16,  17  and  18  show  the  ray  tracing.  The  last  machine  printed 
line  la  the  RADAT  card  for  the  ray  tracing. 

The  makeup  of  the  lens  can  always  be  checked  by  the  input  pi*lnting. 

The  cards  are  all  in  sequence  in  the  deck  even  though  they  appear  to  be 
spread  sqpart  in  the  printing. 

The  output  for  the  energy  distribution  calculation  is  shown  on  pages 
30  and  31  of  the  triplet  example.  The  printing  on  this  sheet  is  self- 
explanatory. 

0.  ALTERING  FIXED  CONSTANTS  IN  TOE  PROGRAM. 

There  are  several  constants  that  the  designer  may  want  to  alter  in  the 
program.  This  may  be  done  easily  by  introducing  a  card  marked  LOAD  in  field 
one.  Field  two  should  contain  a  four  digit  location  number.  Field  three 
is  reserved  for  the  constant.  The  constant  must  be  inserted  with  all  ten 
digits . 

The  IDAD  cards  may  be  inserted  directly  after  a  FORD  card.  Ftor  ex«q>le, 
see  Page  12  of  Sample  3. 
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The  following  constants  are  frequently  dhangedt 

Use  of  Constant  Location 


1)  Maximum  allowable  change  in  a  variable.  7116 

2)  Maximum  number  of  passes.  57^8 

3)  Tolerance  on  aberrations.  The  tolerance  Is 

a  per  cent  of  the  target  value  plus  a  fixed 
increment.  The  per  cent  value  Is  In  9589 

The  fixed  increment  Is  in  9588 

li)  The  constants  a,  b,  c  used  in  clear  aperture 

solve  a  li278 

b  U280 

c  U279 


Programmed 

Value 

4-5020000000 

40000000006 


♦5010000000 

♦U710000000 

U910000000 

5020000000 

5010000000 


H.  OPERATION  OF  THE  MACHINE.  CAJU)  ORIENTATIOM. 

1.  Set  reader  switch  (1,2)  to  A,  B. 

2.  Set  punch  switch  (1,2)  to  B,  B. 

3.  Mount  Ordeals  Tape  on  10. 

U.  Mount  Scratch  Tapes  on  11,  20,  22. 

(20,  22  necessary  only  for  SAM  and  RED.) 

•  5.  Initialize  for  Tape  10. 

6.  Press  Start.  Ehd  of  File  on  Reader,  Staz^  on  Punch. 

7.  Punch  In  time  clock. 

8.  Press  computer  reset  and  start  on  typewriter  console. 

9.  If  typewriter  types  out  CDERR 

1)  Remove  cards  from  stacker. 

2)  Remove  from  hopper. 

3)  Run  out  cards  In  machine. 

U)  Add  to  caxds  In  hopper  and  rerun  by  hitting  start  on  typewilter. 
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10.  In  case  of  overflow  or  uncorrectable  card  errori 

Depress  Computer  Reset  and  Start.  Ibis  then  bypaaaes  the  lane 
and  goes  on  to  the  next  lens.  A  conment  RSTRT  is  printed  out  on 
pxdnt-out  page. 

11.  Upon  completion  of  problem^  the  console  typewriter  prints  out  the 
message. 

Print  tape  11  off  line. 

12.  Use  PEST  output  program  for  lUOl. 

13.  Put  output  tape  on  tape  unit  2. 

II4.  Console  Switch  A  up,  all  others  down. 

15.  Put  X  11  paper  in  lUOl. 

GENERAL  COMMIMTS. 

1.  The  third  and  fifth  order  programs  take  approximately  ^  second  per 
surface. 

2.  Ihe  ray  trace  takes  0.1  seconds  per  surface. 

3.  It  takes  approximately  10  seconds  to  do  the  matrix  algebra  required 
in  a  1^0  X  1^0  matrix. 

1|.  Ihe  size  of  the  matrix  is  equal  to  the  sum  of  the  number  of  VAIQf 
cards  plus  the  number  of  CORR  cards. 

5.  The  use  of  weights  on  the  VARY  cards  is  to  change  their  relative 
Influence  with  respect  to  other  variables.  As  one  Increases  the 
weight  from  0  to  9  on  a  VARY  card,  it  tends  to  decrease  the  change 
In  the  variable.  It  has  been  found  that  for  lenses  with  a  focal 
length  around  10,  it  is  not  necessary  to  Introduce  weights  (l.e., 
weights  of  zero)  on  tho  curvatures  or  thicknesses.  In  long  focal 
length  lenses.  It  Is  usually  advisable  to  put  some  weight  on  the 
curvatures  if  it  Is  desirable  for  them  to  Influence  the  correction. 
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•  Wele^ts  may  also  be  assigned  to  the  HOLDT  cards.  A  large  number  for 
the  weight  places  a  heavy  emphasis  on  the  minimization  feature.  If 
the  weight  Is  too  large.  It  may  make  It  Impossible  for  the  target 
values  to  be  reached.  On  the  other  hand,  If  the  weights  are  not  hl^ 
enough,  the  minimization  may  have  little  affect. 

7»  The  problem  of  introducing  weights  is  very  diffic\ilt  and  our  present 
program  uses  two  methods  for  doing  this.  We  are  not  satisfied  with 
either  method,  so  we  are  not  including  them  In  this  write-up.  Appendix 
three  contains  a  discussion  of  the  weighting  problem,  and  describes 
how  It  is  now  being  done. 

8.  It  Is  advisable  to  limit  the  number  of  iterations  to  six.  Usually 
something  Is  amiss  if  this  is  not  enough.  In  long  systems  it  Is  good 
economy  to  cut  this  to  four.  This  gives  the  designer  a  chance  to 
see  If  eveiything  appears  satisfactory. 

9.  There  is  a  mixed  blessing  written  into  the  program.  This  is  a  maxi¬ 
mum  allowed  change  in  variables.  If  any  change  exceeds  this  value, 
it  is  scaled  to  the  maximum  allowed  change  and  all  other  changes  are 
scaled  by  the  same  amount.  This  feature  prevents  many  blow  ups  when 
solutions  are  far  away  from  the  target,  but  it  also  slows  down  the 
Iteration.  The  value  we  have  inserted  in  the  program  works  well  for 
lenses  with  a  focal  length  around  10.  It  is  much  easier  to  Insert 
this  constant  in  hindsite  than  it  is  in  advance. 

.  There  is  a  limit  on  the  number  of  surfaces  >4ilch  may  be  used.  The 
present  limit  is  30,  including  the  object  and  Image. 
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Q.  H.  Spencer 

Inatltute  of  Optics,  University  of  Rodiester,  Rodiester  20,  New  York 


The  use  of  progr— med  ooaiputlng  nechlnery  for  autonstlo  lens 
correction  requires  e  definite  prescription  according  to  which  a  lens 
qyst—  —y  be  adjudgedt  (a)  acceptable  or  not  acceptable,  or  (b)  improved 
or  not  Ijqproved  over  a  previous  configuration.  Ju^ents  of  the  first 
kind  may  be  made  on  the  b— Is  of  whether  or  not  a  given  set  of  equations 
are  satisfied}  judgments  of  the  second  kind,  on  the  basis  of  whether  or 
not  the  value  of  a  "merit  functl—*'  h—  been  reduced.  A  typical  lens 
design  probl—  will  Involve  both  absolute  requirements,  to  idileh  a  Judgment 
of  the  first  kind  Is  iq)proprlate,  and  relative  requirements  calling  for 
a  Judgment  of  the  second  kind.  This  p^er  describes  a  linearisation  method 
designed  to  accommodate  requirements  of  both  types.  Several  previously 
d— orlbed  Unearlutlon  procedures  are  ^own  to  be  Included  within  the 
framework  of  the  present  method. 
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I.  Introduotlon 


The  dsralopKflnt  of  apeed  progranmad  oomputing  aadiinezy  haa, 
dnring  tha  paat  deoadoj  proaqptad  the  Inraatlgatlon  of  rarloua  aethoda  for 
aatoawting,  or  partially  autoauitlng«  the  lana  daalgn  prooaaa*  The  problaai 
la  ooBg>lloatad,  unfortunately,  by  tha  fact  that  tha  funotlona  idildi  MaaoTa 
qrataa  parfoxaanoa  are  non-Unaar  In  tha  ayatam  rarlablaa.  One  approach  to 
tha  problaa  la  to  approxlaata  tha  non-linaar  funotlona  by  linear  onaa  and 
aolra  the  linear  problaai.  The  Iterative  appUoatlon  of  auoh  a  prooadora  oan 
be  aapactad  to  oonvarga,  under  certain  oondltlona,  to  the  aolutlon  of  the  ^ 
nen-Unear  problaai.  Ihla  iqiproaoh  la  adopted  In  the  preaent  paper. 

In  tha  Method  to  be  deaorlbed,  Lagrange  aultlpllera  are  introdnoad 
aa  a  neana  of  achieving  a  aufflolently  general  framewoilc  to  annnaaaiilatn  both 
abaolute  and  relative  ayaten  requlrenenta.  Abaolute  requlraaienta  are  taken 
hare  to  naan  thoae  reqolrementa  tdiloh  nuat  be  met  exactly  In  order  that  a 
•yataai  be  oonaldered  acceptable.  Relative  requlreaianta  are  thoae  vhloh  may 
be  luaqped  together,  vlth  varloua  welf^tlng  faotora.  Into  a  alngle  tMrlt 
ftuotltm*  ahoae  value  la  to  be  mlnlnlaed.  Sxperlenoe  Indloatea  that  It  la 
dealrable.  If  not  naoaaaazy,  to  allow  requlrenanta  of  both  typea  to  be  lagioeed 
In  the  oourae  of  dealgnlng  a  lena  ayatam. 

It  ahould  be  mentioned  that  the  uae  of  Lagrange  aultlpllera  la  not 
now  to  aatcauitlo  lena  daalgn.  Feder^,  and  Melron  and  Lobenateln^  have  daaorlbod 
modified  ateopeat  deaoent  mathoda  In  vhloh  Lagrange  aultlpllera  are  uaed  to 
detendJie  a  direction  of  change  for  the  eyatam  varlablea  auch  that  certain 
propertlea  of  the  eyatam  remain  unchanged  while  the  value  of  a  merit  function 
la  reduced.  Steapeat  deaoent  nethoda,  however,  poaaeaa  aerlooa  deflolanelea 
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fx^oai  a  praotloal  standpoint  j  the  chief  among  these  being  extremely  alow  con- 
yergenoe.  Ihere  ia  aone  evidence  to  indicate  that  linearisation  procedures, 

Shan  they  converge,  are  likely  to  converge  much  faster  than  descant  nsthoda.^*  ^ 

II.  Statement  of  the  Problem 

One  is  given  an  optical  system  of  J  variables,  pj  (J  ■  1,*..,J). 

Those  variables  may  be  surface  curvatures,  axial  separations  of  surfaces, 
refractive  indices,  or  aspheric  deformation  coefficients.  In  terms  of  those 
variables,  a  number  of  functions,  fj^  (Pj^, •  •  •  ,Pj ) »  "‘■y  defined  idiloh  serve 
to  iseasure  varLoua  diaraoterlstlcs  of  the  system,  both  with  regard  to  its 
physical  structure  and  its  performance.  The  design  problem  is  that  of  finding 
a  slnultanaous  solution,  or  something  approaching  a  simultaneous  solution,  of 
the  set  of  equations 

fj^(Pl, •  ,Pj)  ■  (it  "  (1) 

idiere  the  Sj^  are  constants  representing  the  desired  values  of  the  fj^. 

Aside  from  such  basic  system  properties  ns  focal  length,  magnifi¬ 
cation,  back  focus,  and  total  length,  the  choices  ^Ich  exist  for  the  functions, 
fj^,  are  remarkably  numerous.  Among  them  one  finds  monochiomatlc  and  chromatlo 
aberration  coefficients,  simple  ray  deviations  referred  to  either  the  chief 
ray  or  Oausslan  image  points,  mean  spot  diagram  radii,  weighted  moments  of 
spot  diagram  distributions,  mean  square  wavefront  deformations,  and  energy 
distribution  and  frequency  response  characteristics.  Recent  attempts  to  evolve 
single  "figures  of  mez*lt"  for  Images  of  extended  objects  have  produced  such 
quantities  as  the  relative  structural  content,  fidelity  defect,  and  ooirelatlon 
quality^.  No  attempt  will  be  made  here  to  assess  the  relative  merits  of  these 
various  choices.  It  should  be  pointed  out,  however,  that  considerations  of 
physical  significance  (l.e.,  whether  or  not  a  meaningful  Indlcatlcxi  is  given 
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of  tho  wgr  In  which  a  glYon  object  field  will  be  Inaged)  nuat  be  mpplcMnted 
hgr  ecneidwetioiui  of  ocogmtetionel  effioiencjr  and  llnearitj. 

NoWf  fireqaontly  e  slnultaneoua  solution  will  npt  exist  for  tho 
entire  set  of  eqiuetions.  If  for  exanple,  the  existenoe  of  e  slxnltenooas 

solution  is  extrenelj  unlikely.  In  such  oeseSf  the  adniniBetlon  of 

k-1 


liiere  the  w^  ere  weii^t  feotors^  nay  be  an  eoceptable  altemetiTe  to  an  exeot 
solution.  If  a  slnultaneous  solution  does  in  fact  exist,  it  will  oorrespcnd 
to  the  absolute  artnlimm,^*  0. 

In  general,  there  will  be  oertain  of  the  functions,  for  idhioh 
it  is  essential  than  an  exact  solution  be  obtained,  while  for  the  renainlng 
fonotians  tho  ninisdiatlon  of  ^wlll  eoffloe.  One  thus  dlwldes  the  fj^  into 
two  gnnqw, 


B«(Pl»  • 
and 

Vpi»* 

and  seeks  to  nlnlniBe 


•Pj) 

(n  ■  1, . . .  ,M) 

Pj) 

(n  -  1,..,,M^J) 

M 

^  *2^  ^i^8n  -  *n)^ 

in»l 

subject  to  the  simltineous  solution  of  the  set  of  N  equations, 
where  the  t^  are  the  required  ralues  of  the  funotibns,  h,^. 


(3) 


(U) 


k 


III*  Lin— rlMitl(m  of  Frobl— 


Am  prorlooaly  Indloatod,  tho  funotlona,  and  ara  ganerally 


Bon-Unaar*  Howarar,  in  tha  nalf^borhood  of  a  glran  pointy  (p^«***«p')«  thagr 
■agr  bo  r^praaontad  tgr  tho  oonatant  plua  flrat  ordar  tan—  of  thalr  Tigrlor 
— panalona  about  that  point.  To  thla  approxlnatlon*  than, 

^  (3gn/9Pj)(Pj-Pp  (5) 

and 

J 

h„  -  h^  *»  (^hn/3pj)(pj  -  pp,  (6) 


than  tha  dorlTatlToa  ara  araluatad  at  the  point  (p£,...,p^). 


For  notatlonal  oonTenlanoa,  the  folloulng  deflnltlona  an  now  —dot 
"b  -  8m  -  <>»  ^ 

^  *  Ni  “  ®n  / 

Pj  -  P]  •  ‘Ij  V  (7) 

^8«0pr  So 

3h  ^  b  4 

nj  J 

Ualng  tha  Unoar  approxLmatlona  (5)  and  (6j  for  g,  and  h^^,  Eqa.  (3J 
ahd  (li)  baooaw.  In  tarna  of  tho  deflnltlona  (7), 


J 

^'’nj'lj’^n  (n  -  (Ua) 

i<L 
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TT,  Ae  Method  of  Lagrange  Multipliers 

Hie  probleai  of  finding  an  axtrsana  of  a  function  subject  to  anxillaiy 
constraints  can  be  most  profitably  attacked  bf  the  method  of  Lagrange  multi¬ 
pliers^.  Altfaou^  in  the  present  context  its  ipplioation  mill  be  restricted 
to  the  minimisation  of  a  quadratic  fbnotion  subject  to  linear  oonstraints,  the 
method  is  quite  generally  applicable. 

Let  the  aquations  of  constraint  for  a  systen  of  J  rariables  be 
represented  by 

■  Sjj  (jj  «  1,...,N^J),  (8) 

ahere  the  e^  are  constants.  Ordinarily  such  a  set  of  equations  mill  not 
possess  a  unique  solution.  Rather,  there  will  be  a  continuum  of  points 
satisfying  the  equations.  If  (q^,...,qj)  is  one  such  solution,  then 
(qj^  ♦  dq]^,...,qj  *  dqj)  is  also  a  solution  prorided  the  differentials, 

(di^,...,dqj),  satisfy  the  relations 

J 

■  Z  ^^^J^'^J  *  °  (9) 

J-1 

It  is  required  that  a  solution  point  be  chosen  at  which  a  given 
function, ^(q^,...,qj),  remains  stationary  with  respect  to  differential 

variations  consistent  with  Kqs.  (9).  At  sudi  a  point,  then, 

di-  i(dj/dqjdq.  -  0,  (10) 

*  J-1  J  J 

where  the  differentials  satisfy  Sqs.  (9). 

At  this  Juncture,  the  artifice  of  Lagrange  multipliers  may  be  intro¬ 
duced.  If  Bqs.  (9)  are  satisfied  by  a  set  of  differentials  (which  was  the 
condition  under  idiloh  the  extremum  of  was  defined  above ) ,  then 


6 


n-1  n-1  j-1 

Is  alto  sstlsflsd  bj*  that  sat  of  differentials,  idiere  theXn  arbltraxy 
■nltlplisrs*  Bras  Iq.  (10)  id.ll  be  unchanged  b|y  addition  of  the  sum. 


%• 


PerltoXMdnf  this  addition,  Eq.  (10)  becomes 


(12) 


Now  Iq.  (12)  will  be  satisfied  If 

(5$/^qj^)  °  "  1,...,J).  (13) 

Bqs.  (8)  and  (13)  together  fonn  a  set  of  N  ♦  J  equations  in  the 

N  e  J  unknowns, X  *nd  q.*  and  are  potentially  solvable  for  these  unknows. 
n  j 

Th%  nmltwt  aat  of  Taluosi  qy  will  then  aatlaf/  the  extrenum  condition^ 

Bq«  (10),  In  addltiop  to  satlsf^rlng  the  equations  of  constraint. 


V.  Application  to  the  Present  Problem 
The  present  problen  requires  the  ndnlmlsation  of  the  function, 
given  by  Eq*  (3a)  subject  to  the  constraints  represented  by  Eqs.  (Ua)* 
we  set 


§. 

Thus 


n' 


'V  -Z  %’'y 


3-1 


(111) 


from  uhloh 


(15) 
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Also,  from  Iq.  (3«)« 


Pm  j 


^1 


Iqs.  (13)  thua  booom 


-  •5iS*Sl.-Vi  1 


S  ♦S  "  Ssi  •mk‘V  (k-l,...,j),  (17) 

m-ll^  *'  n-1  m-1 

liioro  V*  hu  bem  replaced 

Iqa.  (17)  and  (Ua)  form  a  aet  of  N  ♦  J  linear  oquatlona  ahioh  nagr 
be  oolTed  bgr  atandard  technlouev  for  the  N  ^  J  unknowia^  qj  and*^* 

Before  proceeding  farther.  It  ahould  be  rerlfled  that  the  extremum 
of ^  guaranteed  br  the  aolutlon  of  Eqa.  (17)  and  (Ua)  la  aotuallj  a  minimum. 
If  { la  expanded  In  a  Tagrlor  aerlee  aboat  the  extremum  point,  the  flrat  order 
tenaa  Tonloh  bf  rlrtne  of  Bq.  (10).  Thua  In  the  nel^borhood  of  the  aitrmaimi, 
I*  llren  by  the  eeoond  order  terma  of  the  Taylor  aerleat 

(xe, 

liiere^,  denotea  the  extromum  ralue.  It  la  oaoumad  that  the  rarlatlcaa, 

A  andiAq|^,  «re  oonalatent  ulth  Iqa.  (Ua). 
from  Iq.  (16), 


00  that  Iq.  (16)  becomeo 
J  J  M 


J»1  k-1  m-1  ^  •'  m-1 1  J-1 


(20) 


8 


IlmSf 

(a) 

rttidb  riww  la  a  aliilana. 

If  tha  aqnalltj  In  Iq.  (21)  holda  for  any  allowad  Tarlatlona,  then  ^ 
la  aald  to  bo  aaBd.«^daflnlta  and  tha  alnljmi  la  not  uniquely  defined.  In  thla 
eaaOf  Iqa.  (17)  aid  (lia)  do  not  poaaeaa  a  unique  aolntion.  The  aabigulty  In  the 
adlatlen  wgr  be  raawredj  howarer,  by  adding  to  the  ana 

J  2 

®  '  (22) 
d-1 


tha  0.  are  aal^t  faotora.  In  addition  to  reacwring  tha  aablgnlty  In  the 
J 

aolatlenf  the  Inclualon  of  the  aua,  S,  afforda  a  neaaure  of  control  over  the 
Inflnanae  of  the  dlfferact  Taxlablea  on  the  aolutlon.  A  large  Talna  for  the 
wallet  faetor  Oj^,  for  aaonqile,  night  be  expected  to  yield  a  relatlrely  aaiall 
aolntion  aalna, 

If  8  la  to  be  Included  In^,  Iqa.  (17)  nnat  be  nodifled  by  adding 

tha  tem 


i  09A^)  -  Cj^. 


(23) 


Iqa.  (17)  are  thoa  replaced  by 


S  Wj  *  Wt  *’■**& ‘S’Si 

n^l  h«l  n«l 


(k  ■  l)...yj).  (2U) 


fl.  Matrix  Ftom  of  t^  Solution 

The  nanlpolatlon  of  8qa.  {Hi  and  (lia)  nay  be  greatly  facilitated  by 
reeaatlag  the  aqpaticna  in  natrlx  fom.  The  neceaaazy  natrlcea  are  tha  ftolloHingt 
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9. 


(25) 


(26) 


(27) 


(28) 


(29) 


10. 


1h«  aatrlx  A  and  ooluan  matrix  d  wLll  almya  appear  ■ultlpllad  bgr 
the  aatrlx  of  val^t  faotora  V.  For  oonvanlanoe  thaaa  produota  aro  daflnad 
to  bo 

M  -  MA  (33) 


f  -  Wd.  (3li) 
AlaOf  In  ovdor  to  oonoorro  opaoo,  a  aatrlx  0  and  a  oolonn  natrlx  g  aro.doflnod 
by 

0  -  mSi  ♦  0  (35) 


g  -  M^f, 


(36) 


lAkoro  tbo  aaporoortpti  t,  donotoo  the  tranepoae. 

With  those  doflnltiona,  Iqa.  (2U)  and  (l^a)  boooao,  roopootlroly, 

Qq  ♦  bO-  g  (37) 


and 


Bq  ■  0 


(38) 


f 
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MiUa  the  eolutlcm  of  Bqa.  (37)  end  (36)  may  be  obtained  throuf^  the 
Inve^Nilon  of  a  alngle  aatrlx  hanrlng  diaianalona  (NxJ)x(NxJ)«  It  is  pexhaps  better 
for  reasons  of  oosputational  aoouraoj  to  arrire  at  the  solution  bgr  inrertlng 
tvD  snaller  natrioes*  Fim  Sq.  (37)* 

q  -  0-^(g  -  _  (39) 

Substituting  this  result  in  Iq.  (38)  and  solving  forT^ 

z)- E"^(BO“^g  -  e),  (J|0) 


lAiere 


E  -  BO 


OnoeVhas  been  detemilnedf  it  may  be  substituted  in  Eq.  (39)  to  determine  q. 

The  matrloes*  0  and  S*  idiioh  required  inversion  have  dimensions  JxJ  and  NxM, 
rsspeotlvely.  Both  of  these  matrices  are  synmetrlc  and  hence  require  fever 
(qieratlons  for  their  inversion  than  do  arbitrary  matrloes. 

Slnoe  the  solution  of  Eqs.  (37)  and  (36)  is  based  on  the  linear 
approximations  expressed  by  Eqs.  ($)  and  (6),  it  vlll  not  in  general  satisfy 
Eqs.  (li)  and  yield  a  true  minimum  of^.  However,  if  a  solution  does  in  fact 
exist  for  Eqs.  (U),  then  Eqs.  (37)  and  (36)  may  be  applied  iteratively  to  arrive 
at  the  required  solution. 

It  is  Inportant  that  the  changes  in  the  variables  produced  at  a  given 
iteration  be  kept  within  reasonable  limits  so  that  the  linear  approximations 
retain  sonM  degree  of  validity.  Otherwise,  a  system  configuration  may  be 
generated  which  is  so  far  removed  from  the  final  solution  that  the  chance  of 
convergence  is  reduced  to  the  vanishing  point.  If  the  e^^  have  been  reduced 
to  aero,  the  weight  factors,  Cj,  may  be  used  to  limit  the  changes  produced  at 

a  given  iteration.  If  the  e^  are  not  all  aero,  however,  it  is  liqposslble  to 
confine  the  changes  in  all  of  the  variables  within  arbitrary  limits  by  use  of 
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th«  waiflittf  0^.  In  suoh  oaaea  one  ougr  take 

(J  •  (1*2) 

itiere  k  la  asaicnad  a  lan,  mcn^  rtUxm  to  prtmob  may  ^  tha  pj  Arm  amwadinc 
qpeelflod  Ualta. 


VII.  Choice  of  Welf^t  Faetora 

The  veleht  faetora,  and  Cj,  play  an  liqportant  role  in  detendning 
both  the  natnra  of  the  aolutlon  obtained  and  the  apeed  with  vhlch  the  wathod 

eocirergea  to  the  aolutlon.  It  haa  already  been  indicated  that  the  faetora  o . 

J 

control  the  effectlTeneaa.  of  the  different  varlablea.  Aa  a  factor  e^  la  wade 
larger,  the  correapondlng  variable,  qj^,  la  forced  toward  a  analler  aolutixm 
value.  Thla  control  ahoold  be  of  particular  value  In  preventing  varlablea  fraa 
atraylng  far  beyond  phyalcally  realizable  lijnlta.  The  factors  detezwine 
the  en^aale  placed  upon  the  corresponding  terma  of  ^ .  A  large  value  for  a 
particular  factor,  w^,  will  lead  to  a  aolutlon  for  ihlch  the  k-th  tern  of  ^  la 
likely  to  be  reduced  to  a  greater  extent  than  If  a  lesser  value  of  were 
used.  The  overall  balance  between  the  values  of  the  factors  w^^  and  the  factors 
0 j .  detendnea  whether  the  major  en^hasls  Is  placed  on  the  minimization  of  the 
quadratic  approximation  to  f  or  on  keeping  the  changes  in  the  syatem  variables 
relatively  small  from  Iteration  to  Iteration. 

B\ren  when  all  the  welj^t  faetora  are  equal,  however,  there  may  be  an 
artificial  wel(^ting  which  arises  front  the  fact  that  both  variables  and  per- 
foxmance  functions  may  be  a  variety  of  types  having  widely  differing  ranges. 

For  example,  a  variation  of  10  in  the  value  of  a  hlg^  order  aspheric  defer-' 
matlon  coefficient  mi^t  produce  dbangea  In  a  set  of  perfoxvance  functions  of 
the  sane  order  of  magnitude  aa  those  produced  by  a  variation  of  10  In  the  .  value 
of  a  lens  elasMnt  thickness.  Thus,  with  equal  weights,  Oy  for  these  two 
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TarlaUM,  a  noold  be  obtained  In  idiloh  virtually  all  of  the  "work” 

eae  dOM  bgr  the  aepherlo  ooefflolant*  Sljrilarljrf  If  one  tena  of  ^  had  a  larger 
value  end  larger  derlvatlvea  than  any  of  the  other  teanaa^  the  greatest  Mpbaels 
would  be  placed  on  the  reduotlon  of  that  taiw. 

'  It  la  of  prluMiy  ijigMrtanoe  that  the  effeota  of  artlflolal  wsi^tlng 
of  the  variables  be  counteracted.  Because  of  the  large  dlfferanoea  aaong  the 
variables  uhioh  swy  be  encountered  in  practice  ^  It  Is  possible  for  certain 
variables  to  be  rendered  alinost  totally  ineffective  hy  artificial  wali^tlng* 

If  significant  (diangea  in  these  variables  are  necessary  to  the  final  solution^ 
the  aathod  My  fall  to  converge  or,  at  best>  will  converge  slowly* 

A  possible  way  in  which  to  counteract  the  artificial  weighting  of 
the  varlablea  Is  to  require  that 


vhsre 


(U3) 


This  yields 


M 


N 


m-l  n»l 


(U5) 


In  order  to  retain  the  ability  to  control  the  welgiiting  of  variables,  Oj  nay 
be  replaced  by 


2  -  «  2 


'j 


y  w*l  n"l  / 


(U6) 


Iqpal  values  for  the  new  wslc^t  factors,  c^,  will  have  the  effect  of  making 
all  of  the  variables  about  equally  effective* 


A  nuiMrloal  problem  remaliiB  in  that  the  elements  of  the  matrices  A 
and  B  may  exhibit  an  extremaly  wide  range  of  values  so  that  round-off  errors 
beeoAs  troublesosw.  This  problesi  may  be  alleviated  by  transfo|«4ng  to  a  nev 
set  of  varlAblea  defined  by 


liiere 


Iqs*  (37)  and  (38)  then  beoome 


idiere 


’3 

(U7) 

o  N  1/2 

m 

TJ  q  g 

m 

(50) 

0  ♦  Z 

(51) 

(52) 

The  matrlossfC^  and  7  are  the  diagonal  matrices  i 

/  -1  V 


0<-l 


\  o 

.  e 

O  •«-! 


o 
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It  la  erldant  trm  ttie  abore  that  the  effeot  of  artlflolal  weightiag 
of  the  Tarlablea  (Haappeara  idth  proper  aealing  of  the  Tarlablaa. 

▼HI*  Partioolar  Oaaea  of  Intereat 
Sereral  preriooaly  tdeaerlbed  aethoda  of  antoeatio  laoa  ooRoetiaii 
appear  vlthln  the  fraMmrk  of  the  preaent  foreulaticai*  It  ahoold  be  inatanio- 
tlTe«  therefore,  to  review  them. 

The  Method  of  Least  Squares  used  by  Rosen  and  Eldert'^,  and  by  Nelron^, 
deals  with  the  mlnljalsatlan  of  ^In  the  absence  of  oonstralnta.  Urns  B,  C,  and 
e  beooM  null  natrloea  and  Iq.  (39)  reduces  to 

q  -  0-^  -  (M*M)-^^.  (55) 

This  method  has  been  used  with  seme  success  to  reduce  the  magnitude 
of  spot  diagram  ray  deviations  and  to  minimize  third  order  aberratlcas.  It  la 
neoesaaxy  that  the  number  of  variables  be  less  than  or  equal  to  the  nuaiber  of 
tezms  of  Otherwise  the  solution  will  be  indeteminate.  In  scow  cases 
eicoesslvsly  large  values  of  the  qj  may  be  produced.  The  changes  applied  to 
the  variables  must  then  be  United  in  the  maimer  indicated  Eq.  (Ii2). 

Pedei^  has  pointed  out  that  idiile  the  method  may  converge  relatively 
quickly  to  the  vicinity  of  a  minimum,  the  accuracy  of  the  method  in  that 
vloliiity  may  be  so  IsqMdred  by  nonlinearity  as  to  prevent  further  developsMnt. 

Hopkins  and  MoCarthy^^.  have  described  a  proceduzv  designed  to 
yield  explicit  values  for  a  set  of  system  functions  when  the  number  of  available 
variables  exceeds  or  equals  the  number  of  functions.  (In  their  applioatlon, 
the  system  functions  were  the  third  order  monochrceiatio  and  first  order  diromatlc 
aberrations.)  The  procedure  follows  from  the  present  fomulatlon  by  eliminating 
^and  setting  Cj  ■  1  for  all  j.  Thus  A  and  d  becoaw  null  matrloes,  0  beooams 
the  identity  matrix,  and  Eq.  (1|0)  reduces  to 


16, 


1/  -  -B“^e  -  ($6) 

Eq.  (39)  then  yields 

q  -  B^(BB^)“^e.  (57) 

This  melhod  has  been  used  for  several  years  at  the  Institute  of  Optics 
and  has  proved  to  be  a  powerful  aid  Ih  practical  deslgn^^.  '  The  condition  that 
the  sum  of  the  squares  of  the  be  minimized  tends  to  confine  successive  changes 
in  the  system  configuration  within  it^asonable  limits  When  the  elesients  of  e  are 
not  excessively  large  nor  the  elements  of  B  excessively  small.  In  oases  for 
liilch  large  changes  have  been  generated^  recourse  to  the  limiting  procedure 
represented  by  Eq.  (42)  has  usually  proved  successful  in  preventing  the  process 
from  diverging. 

Recently,  >(ynne^>  ^  has  described  a  method  of  "Successive  Linear 
Approximatlan  at  Maximum  Steps",  'vdildi  is  an  extension  of  the  Method  of  Least 
Squares.  Here  one  Includes  the  sum  S  but  sets  all  of  the  weif^t  factors,  Cj, 
equal.  Thus 

C  -  cl,  (58) 

idiere  I  is  the  identity  matirlx.  Eq.  (39)  then  becomes 

q  -  (M%  ♦  cl)-lM^f.  (^9) 

This  method  provides  control  over  the  magnitude  of  the  qj  at  each 
iteration  thzough  the  choice  of  a  value  for  c.  As  c  is  increased,  the  magnitude 
of  q  (considered  as  a  J  dimensional  vector)  is  decreased  and  the  direction  of 
q  is  shifted  toward  the  direction  of  the  negative  gradient  of  ^at  the  point 
q  ■>  0,  i.e.,  toward  the  direction  of  maximum  decrease  in  The  ability  to 
restrict  the  q^  to  a  region  of  approximate  validity  of  Eqs.  (6),  coupled  with 
the  removal  of  any  indeteminacy  in  the  solution,  assures  convergence.  It  is 
necessaiy  to  choose  c  carefully  if  the  full  potential  of  the  method  is  to  be 
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rvallscdf  howrer.  Too  large  a  value  uUl  overly  reetrlot  the  diangea  In  the 
varlhblee  eo  that  the  speed  of  oonvergance  will  be  unneoessarily  reduoed.  Too 
■Mil  a  talus  asgr  rasnlt  in  lass  than  optlnua  oonrargenoe  beosnse  of  noo-Unaarllgr. 

n.  Conolnslcns 

It  la  not  antloipated  that  the  nethod  described  here  will  be  a  panacea 
for  all  problasw  of  lens  design.  It  Is  expected,  however,  that  the  nethod  will 
protlde  the  lens  designer  with  a  particularly  flexible  neans  of  eoqploltlng  the 
adrantage  of  high  speed  autoamtle  oalenlatlon  offered  by  nodem  oonputlng  aaohlneiy. 
Ihe  power  of  the  nethod  will  depend  largely  on  the  skill  idiioh  the  dealgnar  dls* 
plays  In  eontrolUng  Its  use. 

The  dioloe  of  systen  evaluation  functions  should  dapend  not  only  on 
the  particular  reqolrenents  of  a  given  design  problen  but  also  on  the  stage  of 
develop— it  of  the  design.  At  the  outset  It  my  be  nost  expedient  to  woric  with 
sijg>le  evaluation  functions,  e.g.,  the  third  and  fifth  order  aberrations,  leaving 
the  use  of  nore  definitive  but  oosputatlonally  nore  conplex  evaluation  functions 
to  a  later  stage.  If  reasonable  values  cannot  ^  obtained  for  the  sinple  fonotiona 
it  will  be  neoessaiy  to  nudce  a  najor  alteration  in  the  systen,  either  by  adding 
aspbearios  or  additional  elenents,  or  by  selecting  an  entirely  different  Initial 
oonfiguration.  ]li  such  drcunatanoes  the  additional  tine  required  to  ocnpute 
values,  for  nore  elaborate  functions  from  the  outset  would  be  wasted. 

It  has  been  mntloned  that  the  choice  of  weight  factors  for  the  tenu 
of^detemines  the  nature  of  the  solutlm  obtained.  As  a  design  progresses, 
the  Judgment  of  the  designer  in  modifying  the  wei^t  factors,  so  that  emphasis 
is  shifted  fron  one  oharaoterlstlo  of  the  systen  perfomanoe  to  another, 
bo  influential  in  detendning  the  quality  of  the  final  design. 
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iNTRoroonow  to  oeometric  optical  BaouBtcr  wppoiisi 


I.  Treneltion  from  the  Wave  Optical  to  the  Oeonetrlo  Optical  Heeponee 
Hie  iMve  optloel  Intensity  response  Is  given  by 

Xni,)  dZdl,  (1) 

shers 

A(X.I)  -  |A(X,I)|e“  (t) 

end 

|a(X,I)|  *  ‘Odl.  (3) 

X  represents  the  wavelength  of  lif^t  In  the  Ijuge  speoei  the  radius  of 
the  Osttssian  reference  sphez^i  Ny)>  the  spatial  frequency  ooaiponontsi 
(XyT)«  the  ooordlnatee  of  a  general  point  on  the  reference  sphere,  referred 
to  coordinate  axes  in  the  exit  pupil  plane* 

A(X,I)  is  called  the  ocsqplex  aperture  tranaadaslon  funotion,  or 
sljg>3y  the  aperture  function*  It  specifies  the  enplltude  and  phase  variation, 
over  the  reference  sphere,  of  the  li(^t  originating  at  a  nonodironatlc  object 
source  point*  |  A(X,T)  |  is  usually  assumed  to  have  a  constant  value  of  unity 
within  the  physical  limits  defined  by  the  ead.t  pupil  and  aero  value  outside 
these  limits*  (Oils  assumption  is,  of  course,  invalid  for  systeaui  having 
apodisiag  ^>erturos*) 

0  (X,I)  is  called  the  wave  aberration.  It  is  defined  with  respect 
to  a  oonstant  phase  wavefront  in  the  nel^bortiood  of  the  exit  pupil*  To  be 
definite,  we  may  take  the  wavefront  which  passes  through  the  origin,  (o,o). 


-1- 


-2- 


at  tha  «d.t  pi«^>il.  (pilfl)  !■  tha  distance  fron  a  point  on  the  narafnixt  to 
the  point  (X»T)  on  tho.rotoronoa  ephoro,  noasurod  along  the  nomal  to  tha 
lanrafiaut  (i.a.i  dLong  tha  rajr  ahldi  intarsaots  the  raferanoa  sphara  at  (XfT)). 

Our  purposa  In  this  aaotlon  Is  to  dadooa  tha  form  of  tha  rosponsa 
fonotlon  in  the  short  vavalangth  Unit,  X  —^0.  Ihls  Is  tha  gaonotrio  optioal 
approxlBatlon. 

For  conronlanoa,  ws  first  writs 


i.  Y 

A(X,I)  A*(X-  XRMjj,  T-  XRNy)  •  C  ^ 
liiara,  fron  (2), 

r-  |a(X,I)|  I  A(X-XBIjy  I-A»lyW  • 


<U) 

(5) 
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Noir  Mjr  be  thoui^t  of  ee  an  operator  which  Is  applied  /vw  tiaes  in 

the  nth  tens  of  the  aeries.  We  know  that 

•  -R(W, 'jS  > 


SO  that  the  operator  la 

Henoe  we  nay  write  the  series  (9).  as 

-  #(».n -x^Kir  ,.i 


(10) 


♦l»l 


(11) 


In  the  preceding  equatloni  we  have  made  use  of  the  fact  that  (K  )  ■  (X,T) 
when  X  "  o  *0 


__  A  <  \ 

'  VitX'l'i*') 

\x 

X  V 

khare  9  ( «* ,  ^ )  is  any  function  of  «<  and  ^  . 


Subatltutlng  (11)  In  (6)t 


•»||a 


f 


Ihui, 


I.W  >  ^  ^vR  (  ^  ^  )  • 

X-fco 


(12) 


(U) 


It  la  aaauBMd  that  nona  of  tha  darlTatlvaa  baooaM  Inflnlta.  lhara  ara 
unuaual  oaaaa  In  which  this  aaaumptlon  la  violated  (a»g.f  whan  a  aava  nonMl 
la  tangent  to  tha  rafaranoa  apfiara),  but  the  raault  (13)  reaialna  ganaralljr 
valid. 


Va  ara  nov  able  to  write  tha  geonetrlc  optical  reaponaa  ualng 
tha  raaulta  (13)  md  (7)i 


\-^9 


O*  CkO 


“■'U 


j\jy .  (1^*) 


|.(i,i)|  .{j:  ssWi's:.] 

then  tha  nomallalng  factor,  K,  la  aiinplj  the  area,  a,  of  the  pupil,  and 
tha  Infinite  Unite  in  (lU)  nagr  be  replaced  by  the  Unite  defined  by  the 
pupil.  Thus 
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II. 


Repreaentation  In  Terns  of  linage  Plane  Ray  Coordinates 
Our  purpose  in  this  section  Is  to  relate  the  wave  aberration,  ^  , 
to  the  Image  plane  coordinates  of  the  rays  associated  with  the  given  wave- 
front.  Ihe  result  will  be  used  In  Eq.  (15)  to  obtain  an  esqpresslon  for 
the  geonetrlc  optical  response  which  depends  only  upon  the  Image  plane  ray 
cooidlnates.  The  advantage  of  this  representation  comes  from  the  faot  that 
It  Is  an  easier  task  to  calculate  accurate  ray  coordinates  than  to  oaloulate 
wave  aberrations. 

Fig.  1  shows  the  exit  pupil  plane,  (X,Y),  the  Oaueeian  reference 
sphere,  3,  the  given  wavefront,  W,  and  the  image  plane,  Ihs  Oaussian 

ireferenoe  sphere  is  assumed  to  pass  through  the  origin  at  the  exit  pupil 
and  to  have  Its  center  of  curvature  at  the  Oaussian  image  point, 

In  the  image  plane.  The  Image  plane  is  assumed  to  be  separited  from  the 
exit  pupil  by  a  distance  D. 

Now  the  distance  frcmi  a  point  (X,I,Z)  on  the  wavefront  to  the 
Oaussian  Image  point  is 

R'  -  [(X-  Xk)^  ♦  (Y- 

If  we  define  A  to  be  the  distance  from  the  point  (X,I,Z)  on  the  wavefront 
to  the  reference  sphere,  measured  along  the  normal  to  the  reference  sphere, 

we  nay  write 

R'  -  R  ♦  ^  (17) 

where  R  Is  the  radius  of  the  reference  sphere  and  Is  given  by 

«■  r*!.**^*”']  *  (i") 

Frm  (16)-(16),  we  find 

2RA  ♦  ♦  Y^  ♦  Z^-  2(«jj  ♦  l^k  * 


(19) 


FIOURE  I. 


MltHT 


ASSOCIATED  WITH  POINT  OBJECT,  AS  IT 
ear  at  THE  EXIT  PUPIL. 
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NdfW  tgrplealljr. 


r 


io“^. 


•o  to  a  yxy  good  approximation  we  max  the  term  in  (19)*  The 
roault  max  ^  taken  aa  the  equation  of  the  Mnrofx^t  and  It  mill  be  oon- 
▼oniant  to  write  It  in  the  fom 


(20) 


F(J[,I,2)  ■  ♦  ZD-  ♦  Z^)  ♦  -  0. 

Letting  (X]j>T|^)  repreaent  the  Image  plane  ooordinatea  of  a  ngr 
from  (Z(T|Z)  on  the  waTofront^  the  ayienetrlo  equatlcme  for  the  rax 


lit-  X 


'iF 

Tt 


D-  Z 


Svaluatlng  the  derlvatlveai  we  find 


Xk-X 


W 


Xk-  X  ♦ 


•  1 


\.X<- 


(21) 


(22) 


(23) 


(^  ie  a  function  of  X  and  I  onlxi  once  It  la  given,  the  Z  coordinate  on 
the  wavefront  la  determined  bx  (21)).  From  (23)  we  Immediatelx  obtain  the 
relatione 

yXk  ^  Xk-Xk* 


(21») 


TO  a  good  approximation  we  max  replace  4  bx  0,  the  wave  aberration 


Introdooed  In  eeotlon  I.  The  order  of  magnitude  relation  aaeoolated  with 


7 


thla  approxliMtlon  la 


Thus 


(25) 


(26) 


Theae  relatione  may  be  aubatltuted  In  (1^)  to  obtain  the 
required  e3q;>reaelon  for  the  reapcmaet 


(27) 


One  of  the  moat  direct  ways  of  evaluating  the  above  Integral 
la  to  aum  up  spot  dltigram  data.  One  traces  a  large  number  of  rays  throu(^ 
the  optical  ayatemj  arranging  them  such  that  they  are  uniformly  distributed 
over  the  entrance  pupil  (and  hence  nearly  uniformly  distributed  over  the 
exit  pupil).  The  response  Integral  (27)  Is  then  replaced  by  the  sunnatlon 
formula 

r,(N.  ^  irfA-.C rx.v .+  w,(»Y.v  ] 

J  ’  '  JSI  *  •' 


.  J. 

T  J  ^  J 

J*' 


(28) 


where  the  subscript  j  Identifies  the  jth  rayi  and  T  represents  the  total 


number  of  rays. 
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Arereged  ClreulT  Beepooae 

Snppoie  ve  are  oonoanted  with  Ijuglng  a  almaoldal  target  haring 
apatlal  fmqaanojr*  N*  orlantad  at  an  angla  0  vlth  raapeot  to  tha  Z  axis  In 
tha  dbjaot  plana  (saa  Fig.  2).  lha  spatial  fraqoaniqr  ooagwnants  than  ara 


Nx  -  N  sine 
My  -  N  oose 


V 


(29) 


If  ws  tnnsfoin  tha  araqr  darlatlons  (  Sx^,  to  polar  fontf 


fv. 


(30) 


aa  Mgr  arlta  (27)  as  followai 


ir  iM  X(, 

«  jxjy 


exit  pupil 


■ 


4XJY 


(31) 


exit  p\ipll 

lha  ararage  of  (31)  over  all  target  orientations,  e,  constitutes 
a  useful  criterion  of  system  performance.  We  shall  denote  this  average 
hy  f^(N)  and  call  It  the  averaged  circular  response.  Using  (31), 


®  e*'if  I 


’  ^  \  JV  jy  JY 


txif  I  o 

,ir/t 


•hSS  V  {tvM  f^n  i'^t)  JxJY 

I  0 


(32) 
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Now 


^  \  c*s  (X  i?*»  v') 


(33) 


liiMra  Jg(x)  IM  MTo  ordar  Baaaal  i^otlon.  Our  rasponaa  la  Ihua 


f  (M)  - 

axit  pupil  _ _ _ I 

If  apot  dlagraai  data  la  avallabla,  aa  magr  uaa  tha  auHMtlon 


fanaola  ^ 

J*' 


(3Ji) 


(35) 


ahara  T  la  tha  tmbur  of  ragra  and  (ftk) 


i 


la  the  diatanoa  trxm  tha  Qauaalan 


iMfa  point  to  tha  point  of  intaraaotlon  of  the  jth  ray  with  tha  laaga  plana* 
It  la  ooMmon  praotloa  to  praaant  apot  diagram  data  In  the  form 
of  radial  energy  dlatrlbutlon  curvaa .  If  wa  vara  to  canter  a  aufflolantly 
aaall  olroular  apartura  at  tha  Qauaalan  Image  point,  a  certain  number  of 
raya  would  fall  outalda  and  be  blocked.  If  the  radlua  of  tha  aperture  vara 
allowed  to  Inoraaaa,  wo  would  find  more  and  more  raya  paaalng  through  until 
erontually  all  of  tha  raya  wore  admitted.  Tha  radial  energy  dlatrlbutlon 
ourra  la  a  plot  of  tha  fraction  of  tha  total  number  of  raya  paaalng  throu^ 
auoh  an  apartura  aa  a  function  of  tha  apartura  radlua.  laoh  ray  repTCaonta 
the  Uf^t  anargy  paaalng  throuf^  an  alemental  area  of  the  apartura  ao  that 
a  fraction  of  tha  total  number  of  raya  may  be  Interpreted  aa  a  fraction  of 
tha  total  light  anargy  paaalng  throu^  tha  optical  ayatam. 
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FroM  •  radial  enargjr  dlatrlbutlon  ourva  we  may  aalaot  a  aat  of 
polnta  avMh  aa  tha  foUowliigt 


£ 


I 

I 

I 


1 

I 

I 

I 


C.  «  'li 

A  T* 


I 

» 

I 


H 


liiara  la  tha  nunbar  of  raya  paaalng  through  an  apartura  of  radlua 
(f  f  eantarad  at  tha  Oauaalan  laaga  point,  and  T  la  tha  total  nuabar 

of  ragra  In  tha  apot  dlagrant  lha  polnta  art  ohoaan  auch  that 
(aaa  flf.  3). 

If  tha  polnta  art  not  too  vldaly  aaparatad,  wa  nay  aay  to  a  fair 
agiprailMtlon  that  tha  apot  dlagran  oonalata  of 


ray  polnta  hairing  radlua 


a  »  »  « 


I 

I 

I 


I 

I 

I 


I 

m,  n  n  •*  ii 


•  — 

I 

I 

t 

I 

I 

I 

I 

1. 


FIGURE  2 


•INUSOIOAL  TARGET  IN  THE  OBJECT  PLANE  . 
FREGUENCY  COMPONENTS:  NK*NtinB,  N,«NeotB 


RELATIVE 


TYPICAL  RADIAL  ENERGY  DISTRIBUTION  CURVE  SHOWING  VALUES 
USED  FOR  APPROXIMATE  AVERAGED  CIRCULAR  RESPONSE  CALCULATIONS. 
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nms,  (35) 

Mgr  ba  aritten 

'  *»» 

T. 

(36) 

abon 

4a,.  • 

(37) 

and 

(38) 

Kotlng  that 

— •  »  E.  -  E. 

T  * 

(39) 

wo  roMrlta  (36)1 

f(w^»  ^  ^.EiT.tiTWCXf.V  ] 

(Uo) 

ahoro 

'(la) 

1h«  fonnila  (1^0)  requlrea  mudi  lees  computljig  tine  than  does  (35) •  Typioallx^ 
a  ralua  of  20  for  M  la  eufflolent  to  gire  acceptable  aceuraojr  for  design 
oorkt 


Sorias  Kapanslon  of  the  Averaged  Circular  Response  — • 
Moments  of  the  Spot  l?iagrain  Distribution 


A  typical  tenn  of  the  series  (35)  may  bo  expanded  to  yield 


m 


>12 


Dma 


tdoa  iZ.'L 


).,  r..  ('•'■)*• 


^v— 


K  «  N 


Xm0 


(10) 


iAmt* 


and 


m 


m 


Iha  •(i^  naj  ba  racognlaad  aa  tha  avan  aoaianta  of  tba  radial  spot 
dlagnn  dlatrlbutlon.  Raoant  applloatlona  of  automatic  oomputlng  taohnlquaa 
to  tha  daalgn  of  optical  aTStama  hava  baen  baaad  on  aa  a  maasura  of 
tjwtm  parfonuuioa*  Va  aaa  that  •(  ^  datarminaa  tha  raaponaa  in  tha  low 
fraqoanox  ration,  and  that  mlniniaatlon  of  o<  ^  for  a  giran  Imaga  will  naxl- 
■iaa  tha  low  fratiaanojr  raaponaa.  Dila  gain,  howorar,  auqr  ba  at  tha  aaorifioo 
of  good  hifh  frwqamtey  raaponaa.  Iha  following  araapla  iUuatrataa  thla 
point. 

Soppoaa  wo  hawa  an  Inaga  dtaraotariaad  bgr  tha  oqination 

<w 

aharo  ^  la  tha  radial  oadt  pupil  ooordlnatai 

^  -  (X*  ♦  I*)  ^  . 


(U7) 
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lh«  lirat  tani  of  (1|6)  rapraaonta  dafaot  of  fooua,  whlla  tho  aebond  tax* 
rapraaanta  ondaroorraetod  aavonth  order  apharloal  abarratlaii*  Vfa  ahall 
take  the  ooefflolenti  0|  of  the  first  tom  of  (1{6)  as  a  variable  which 
alloaa  ua  to  alter  the  image  distribution.  Variation  of  o  oorraaponda  to 
a  movanant  of  tho  Image  plane  along  the  Z  axis. 

Fig.  (1()  shows  the  response  ^'(N)  for  several  values  of  e, 
aaauming  an  aperture  defined  by 

^  <  1.  (liB) 

The  curve  labelled  A  Is  obtained  idien  c  ■  o.  The  response  is  seen  to  drop 
fairly  rapidly  in  the  low  frequency  region  to  a  value  of  about  0.6  and  than 
to  level  out|  extending  well  into  the  hi^^er  frequency  region  without 
dropping  below  O.U*  Fine  detail  in  the  object  Will  thus  be  preserved  in 
the  Image  I  although  at  reduced  contrast. 

Tho  curve  labelled  6  is  obtained  when  c  ■>  .OU.  This  is  the  value 
of  c  fbr  which  e^^  is  minimum.  The  response  in  the  low  frequency  region 
has  been  oonsidarably  improved.  In  the  hl^er  frequency  region,  however, 
tte  responsa  now  drops  quite  rapidly  to  aero  and  becomes  sli^tly  negative, 
indiaatlng  "spurious  resolution". 

The  curve  labelled  C  corresponds  to  an  intermediate  image  plane 
position,  0  ■  .022.  The  improvement  over  curve  A  extends  well  into  the 
mid  frequency  region.  Curve  C  is  a  decided  improvement  over  curve  B,  the 
slight  raduotion  in  low  frequency  response  being  more  than  compensated  for 
hy  the  large  elevation  in  the  higher  frequency  response.  (In  special  appll- 
oations,  of  oourso,  curve  B  may  be  preferred  for  specific  reasons  detemined 
by  the  applioation. ) 
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